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ABSTRACT 


The  Mei-Yu  front  is  a  quasi-stationary  belt  of  heavy 
precipitation  imbedded  in  the  summer  monsoon  trough  over 
East  Asia.  This  work  used  a  two  dimensional  frontogenesis 
model  to  study  the  Mei-Yu  front,  especially  its  associated 
low-level  jet  (LLJ)  which  often  occurs  during  intense 
convection.  Several  numerical  experiments  were  carried  out 
to  simulate  the  quasi-steady  frontal  structures  in  different 
environments.  The  results  resemble  many  observed  features, 
such  as  the  stronger  temperature  and  moisture  gradients  in 
the  midlatitude  fronts  and  the  stronger  horizontal  shear  and 
LLJ  in  the  subtropical  fronts.  The  LLJ  is  most  conspicuous 
in  the  subtropical  simulation  when  high  humidity  and  surface 
fluxes  are  included.  The  results  further  suggest  that  the 
LLJ  is  developed  through  the  Coriolis  torque  that  is  exerted 
by  the  low-level  poleward  branch  of  a  "reversed  Hadley  cell" 
equatorward  of  the  front.  This  thermally  direct  cell  is 
different  from  the  normal  cross-frontal  secondary 
circulation.  Its  development  depends  on  the  occurrence  of 
intense  convection. 
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I.   INTRODUCTION 

A  review  of  recent  observational  studies  of  the  main 
phenomena  of  this  study,  the  Mei-Yu,  and  the  associated 
disturbances  will  be  presented  first.  There  follows  a 
review  of  the  pertinent  theoretical  studies.  A  description 
of  the  objectives  of  this  numerical  study  of  the  Mei-Yu 
front  will  then  be  given. 

A.   OBSERVATIONAL  STUDIES 
1.  Mei-Yu 

Over  eastern  Asia,  the  transition  period  from  the 
northeast  winter  monsoon  to  the  southwest  summer  monsoon 
possesses  a  distinctive  climatological  phenomenon  which  is 
important  for  regional  weather  forecasting.  The  seasonal 
rainfall  reaches  a  maximum  during  this  period  (Fig.  1.1). 
The  rainfall  during  this  period  is  called  "Mei-Yu"  (plum 
rain)  in  China  and  "Baiu"  in  Japan.  During  this  period,  the 
rainfall  may  be  continuous  or  intermittent  for  several  days 
to  a  few  weeks  and  is  characterized  by  frequent  rain-showers 
and  thunderstorms,  with  rainfall  rates  of  up  to  a  few 
hundred  millimeters  per  day  (e.g.,  Chen,  1983;  Ninomiya, 
1980).  The  Mei-Yu  season  first  appears  in  southern  China  in 
May,  and  influences  the  weather  over  southern  China  and 
Taiwan  from  mid-May  to  mid-June  (Chen  and  Chi,   1980).    It 
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then  migrates  over  the  Yangtze  Valley  and  Japan  from  mid- 
June  to  mid-July. 

Synoptical ly ,  the  Mei-Yu  rainfall  is  associated  with 
the  repeated  occurrence  of  a  front  that  develops  in  the 
mid] at i tudes  and  slowly  moves  southeastward  to  a  quasi- 
stationary  position  extending  from  southern  Japan  to 
southern  China.  Satellite  pictures  reveal  a  nearly 
continuous  cloud  band  accompanying  this  front,  which  is 
often  called  a  Mei-Yu  front  or  a  Baiu  front.  The  frontal 
position  coincides  with  the  early  summer  monsoon  trough  in 
eastern  Asia  and  its  existence  may  be  viewed  as  a 
manifestation  of  the  trough.  The  beginning  of  the  rainy 
season  also  signifies  the  onset  of  the  northern  summer 
monsoon  in  this  part  of  the  world.  The  Mei-Yu  often 
undergoes  sequential  jumps  between  May  and  July  (Fig.  1.2; 
Tao  and  Chen,  1987)  depending  upon  the  strength  and  the 
position  of  the  North  Pacific  Ocean  subtropical  high  (Tao 
and  Ding,  1981).  As  the  Mei-Yu  onset  occurs  in  the  Yangtze 
River  region  in  June,  the  Indian  monsoon  is  also  being 
established  and  the  westerly  jet  stream  begins  to  retreat  to 
the  north  of  the  Tibetan  Plateau.  At  the  same  time,  the 
subtropical  high  in  the  western  Pacific  area  moves  northward 
to  20-25°N.  The  Mei-Yu  season  lasts  through  mid-July,  when 
the  subtropical  high  moves  northward  to  30°N.  At  that  time, 
the  rain  belt  is  to  the  north  of  the  Yellow  River,  which 
signifies  the  beginning  of  heavy  rainfall  period  in  northern 
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China  and  the  ending  of  the  Mei-Yu  in  subtropical  China  (Tao 
and  Chen,  1987) . 

2.  Mei-Yu  Front 

The  climatology  of  the  Mei-Yu  front  over  southern 
China  and  the  adjacent  areas  was  examined  by  Chen  and  Chi 
(1980)  for  the  periods  of  15  May-15  June,  1968-1977.  They 
found  that  the  Mei-Yu  f rontogenes is  takes  place  in  the 
subtropical  part  of  southern  China  to  the  south  of  ~35°N. 
All  the  Mei-Yu  fronts  were  formed  in  the  area  of  20°-35°N 
and  100°-130°E  (Fig.  1.3).  The  lifetime  of  individual  cases 
of  Mei-Yu  fronts  ranges  from  3  days  to  22  days  with  an 
average  of  8  days.  The  kinematic  background  for  the 
f rontogenes is  is  the  deformation  field  between  the 
subtropical  high  over  the  Pacific  and  the  migratory  high 
over  central  China.  They  suggested  that  the  sea-surface 
temperature  contrast  in  the  nearby  South  and  East  China  Seas 
may  also  be  a  favorable  factor  for  the  f rontogenes is . 

In  contrast  to  the  polar  front,  the  temperature 
contrast  across  the  Mei-Yu  front  is  typically  small,  but 
there  is  always  a  marked  moisture  gradient  near  the  frontal 
zone  (Akiyama,  1973a;  Chen  and  Tsay,  1978;  Matsumoto  et-al., 
1970,  1971;  Tao,  1980).  An  analysis  of  the  kinematic 
structure  of  the  Mei-Yu  frontal  system  over  Southeast  Asia 
was  carried  out  by  Chen  (1978).  The  Mei-Yu  frontal  system 
possesses  a  remarkable  shear  line  which  roughly  coincides 
with  a   trough  at   850  and  700  mb  ,   although  the  associated 
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geopotential  gradient  is  rather  weak  in  the  immediate 
vicinity  of  trough.  The  shear  line  approximately  coincides 
with  the  maximum  north-south  gradients  of  mixing  ratio  and 
equivalent  potential  temperature.  To  the  south  of  the  shear 
line,  a  low  level  jet  (LLJ)  that  is  characterized  by  maximum 
mixing  ratio  and  high  equivalent  potential  temperature  is 
found  near  700  mb  (Fig.  1.4).  .  Strong  cyclonic  vorticity, 
horizontal  convergence  and  upward  motion  are  observed  along 
and  to  the  south  of  the  shear  line.  The  cloud  band 
approximately  coincides  with  the  area  of  maximum  cyclonic 
vorticity,  horizontal  convergence  and  upward  motion. 
Intense  convective  clouds  in  the  southern  portion  of  the 
cloud  band  tend  to  be  associated  with  the  maximum  values  of 
these  kinematic  parameters  as  well  with  relative  humidities 
greater  than  80%  in  the  lower-middle  troposphere.  Clear 
skies  over  the  East  China  Sea  and  its  vicinity  are 
accompanied  by  downward  motion  of  dry  air  (relative  humidity 
less  than  30%  at  850  mb ,  less  than  10%  at  700  and  500  mb ) . 

Chen  and  Chang  (1980)  studied  the  structure  and 
calculated  a  vorticity  budget  for  a  Mei-Yu  system  during  the 
period  10-15  June  1975.  During  the  mature  and  decaying 
stages  of  the  monsoon  trough,  the  structure  of  the  eastern 
section  (southern  Japan)  and  the  central  section  (East  China 
Sea)  is  different  from  that  of  the  western  section 
(southeastern  China).  In  the  eastern  and  central  section, 
the   structure   resembles   a   typical   midlatitude   baroclinic 
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front  with  a  strong  vertical  tilt  toward  an  upper-level  cold 
core  low  and  a  strong"  horizontal  temperature  gradient.  On 
the  other  hand,  the  western  section  resembles  a  tropical 
disturbance  with  an  equivalent  barotropic,  warm  core 
structure,  a  relatively  weak  horizontal  temperature  gradient 
and  a  rather  strong  horizontal  wind  shear  in  the  lower 
troposphere.  In  the  mature  stage,  the  upward  motion  reaches 
maximum  values  at  500  mb .  It  is  located  over  the  850  mb 
monsoon  trough  in  the  western  (W)  section,  4°  latitude  south 
of  the  trough  in  the  central  (C)  and  3°  latitude  southeast 
of  the  trough  in  the  eastern  (E)  section.  On  the  north  side 
of  the  trough,  downward  motion  prevails  with  a  maximum  at 
400  mb  in  sections  C  and  E,  and  at  700  mb  in  section  W. 
This  well-organized  vertical  velocity  couplet  is  indicative 
of  a  well-developed  secondary  circulation  normal  to  the 
trough.  In  the  decaying  stage,  the  circulation  becomes 
substantially  less  vigorous  in  sections  W  and  C,  but  it 
maintains  intensity  in  section  E.  Tsay  and  Chen  (1980) 
analyzed  the  processes  which  force  the  vertical  motion  in 
sections  W,  C  and  E.  They  found  the  eastern  portion 
partially  maintains  the  characteristics  of  a  midlatitude 
system,  while  the  western  portion  behaves  similar  to  a 
tropical  system. 

The  Mei-Yu  front  can  be  several  thousand  kilometers 
long,  and  therefore  it  is  a  large-scale  monsoonal 
phenomenon.   However,  the  cross-front  scale  is  considerably 
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smaller,  which  is  within  the  mesoscale  range.  Among  the 
many  mesoscale  features  that  are  observed  near  the  Mei-Yu 
frontal  zone,  the  most  important  are  the  organized  mesoscale 
convective  systems  embedded  within  the  broad  Mei-Yu  cloud 
band,  which  lies  on  the  warm  side  of  the  front  (Fig.  1.5). 
These  convective  systems  tend  to  move  along  the  front  from 
west  to  east  and  produce  local. ly  heavy  rainfall  rates  of 
>  100  mm  day"1  (Akiyama,  1974,  1978;  Chen,  1977;  Ninomiya 
and  Akiyama,  1971,  1972,  1973,  1974). 
3 .   Low  Level  Jet  (LLJ) 

To  the  south  of  the  Mei-Yu  (or  Baiu)  front,  there  is 
often  a  strong  ageostrophic  low  level  jet  (LLJ)  which  is 
usually  observed  between  1.5  and  3  km  (850-700  mb )  and  is 
strongly  correlated  to  the  heavy  rainfall  (Akiyama,  1973b; 
Chen,  1977,  1986;  Matsumoto,  1972;  Matsumoto,  et  al.,  1971; 
Ninomiya  and  Akiyama,  1974;  Tao  and  Chen,  1987).  Akiyama 
(1973a)  observed  that  the  heavy  rainfalls  in  the  Japan  Baiu 
season  occur  most  frequently  to  the  north  of  the  low  level 
jet  stream  and  seldom  occur  in  areas  more  than  200  km  south 
of  the  jet  axis.  Matsumoto  (1973)  found  that  the  peak  of 
heavy  rainfall  precedes  by  6  h  the  peak  of  the  low  level  jet 
speed.  He  also  found  that  the  peak  of  heavy  rainfall  is  in 
phase  with  the  maximum  acceleration  of  the  LLJ. 

In  the  Great  Plains  and  Midwest  areas  of  the  United 
States,  the  low  level  jets  are  a  cl iraatological  feature  with 
the  maximum  frequency  during  the  spring  and  summer  months 
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(Bonner,  1968).  The  low  level  jets  in  this  'area  are 
basically  a  boundary  layer  phenomenon  that  are  characterized 
by  a  diurnal  oscillation.  The  jet  reaches  maximum  intensity 
by  early  morning,  and  is  associated  with  a  nocturnal 
temperature  inversion  (e.g.,  Blackadar,  1957;  Bonner,  1968; 
Holton,  1967;  Lettau,  1967;  Wexler,  1961).  Different 
theories  (discussed  later  in  section  B)  have  been  proposed 
to  explain  the  generation  of  the  low  level  jet  and  the 
seasonal,  temporal  and  geographic  variations.  In  the  British 
Isles,  a  low  level  jet  in  conjunction  with  midlatitude 
cyclones  was  studied  by  Browning  and  Harrold  (1970)  using 
Doppler  radar  observations.  Browning  and  Pardoe  (1973) 
found  that  for  six  cases  of  ana-cold  fronts  occurring  in  the 
colder  months  in  British  Isles,  the  low  level  jet  had  a 
maximum  velocity  in  the  layer  900  to  850  mb  of  25  to  30  m  s1 
and  occurred  during  night  or  day. 

4 .  Mesoscale  convective  systems  (MCS's) 

Rainfall  characteristics  over  Taiwan  and  vicinity 
were  studied  by  Chen  (1977)  and  Chen  and  Tsay  (1978)  using 
dense  hourly  rainfall  reports,  radar  echoes  and  satellite 
pictures.  The  low  level  southwester lies ,  especially  the 
LLJ,  tends  to  create  potential  instability  to  the  warm  side 
of  Mei-Yu  front.  Continued  large-scale  ascent  then  leads  to 
the  release  of  potential  instability  through  organized 
mesoscale  convective  systems.  The  mesoscale  convective 
systems  embedded  within  the  broad  Mei-Yu  cloud  band  can  be 
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of  several  different  sizes  (meso  -a,  -$,  or  -y) •  A  case 
study  by  Chen  and  Tsay  (1978)  revealed  the  existence  of 
these  mesosys  terns  in  Taiwan  with  wavelengths  of  60,  85  and 
170  km  and  the  corresponding  periods  of  2.6,  3.5  and  7.1  h. 
A  spectrum  analysis  (Tao,  1980)  of  radar  echoes  in  the 
rainbands  along  the  Mei-Yu  front  over  the  Chinese  mainland 
also  showed  the  existence  of  organized  mesoscale  convective 
systems  with  wavelengths  of  150,  84  and  50  km.  Akiyama 
(1978)  found  that  the  echoes  within  the  rainband  in  the  Baiu 
season  in  Japan  are  organized  into  mesoscale  clusters  with 
spatial  intervals  on  the  order  of  100-150  km. 

Simulations  of  these  mesoscale  convective  systems 
require  a  three-dimensional  model.  In  our  two-dimensional 
simulations,  we  will  not  study  these  phenomena. 

B.   THEORETICAL  STUDIES 
1  .   Front 

Quasi-geost rophic  frontal  theories  (Stone,  1966; 
Williams,  1968;  Williams  and  Plotkin,  1968)  showed  that 
large-scale  deformation  fields  can  lead  to  surface 
f rontogenesis  and  result  in  strong  horizontal  temperature 
gradients  that  are  characteristic  of  fronts.  However,  the 
structure  of  the  front  and  the  associated  secondary 
circulation  do  not  tilt  with  height.  It  also  lacks  cyclonic 
wind  shear  at  the  front  and  has  several  other  unrealistic 
features,  including  the  infinite  time  required  for  the  front 
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to  form.  Later  studies  using  primitive  equations  (Williams, 
1972)  and  the  semi-geostrophic  equations  (Hoskins  and 
Bretherton,  1972)  removed  these  unrealistic  features  by 
including  the  ageostrophic  advections  of  temperature  and 
momentum.  In  particular,  Hoskins  and  Bretherton  (1972)  used 
the  set  of  semi-geostrophic  equations  to  show  that  the 
nonlinear  effects  of  the  ageostrophic  advections  can  be 
represented  by  the  quasi-geostrophic  solutions,  if  these 
solutions  are  written  in  the  geostrophic  coordinates  and 
then  transformed  to  the  physical  space.  Later,  the  semi- 
geostrophic  theory  was  expanded  to  three  dimensions  by 
Hoskins  (1975)  and  Hoskins  and  Draghici  (1977). 

Fractional  and  diabatic  effects  also  have  been 
studied  by  a  number  of  investigators.  Williams  (1974)  found 
that  vertical  mixing  can  determine  the  slope  of  a  front, 
and,  together  with  horizontal  mixing,  can  determine  the 
cross-width  of  the  frontal  zone.  Keyser  and  Anthes  (1982) 
showed  that  the  proper  simulation  of  the  concentrated  upward 
motion  in  the  surface  frontal  zone  depends  on  the  planetary 
boundary  layer  parameterization.  Williams  et  al.  (1981), 
Ross  and  Orlanski  (1978),  and  Hsie  and  Anthes  (1984)  found 
chat  latent  heat  release  intensifies  both  the  upper  level 
westerlies  and  the  low  level  easterlies  poleward  of  the 
surface  front,  and  the  ageostrophic  secondary  circulation 
across  the  front. 
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2 .  Symmetric  Instability 

Within  extratropical  cyclones,  banded  structures  in 
clouds  and  precipitation  have  been  observed  both  in 
conditionally  stable  and  unstable  atmospheres  (Hobbs,  1978). 
Bennetts  and  Hoskins  (1979)  suggested  that  these  bands  are  a 
result  of  conditional  symmetric  instability.  Emanuel  (1979) 
demonstrated  that  circulations  resulting  from  this 
instability  are  fundamentally  mesoscale  in  character. 
Recent  observations  indicate  that  conditional  symmetric 
instability  is  responsible  for  some  precipitation  bands 
within  extratropical  cyclones  (Bennetts  and  Sharp,  1982; 
Moore,  1985;  Parsons  and  Hobbs,  1983;  Sanders  and  Bosart, 
1985) . 

Bennetts  and  Hoskins  (1979)  and  Emanuel  (1982) 
assessed  the  degree  of  conditional  symmetric  instability 
based  on  infinitesimal  displacements  in  saturated 
atmospheres.  Emanuel  (1983a, b)  considered  parcel 
instability  which  leads  to  a  measure  of  the  stability  of  the 
moist  baroclinic  atmosphere  to  finite  slantwise  reversible 
displacements  of  a  two-dimensional  air  parcel  or  air  "tube". 
Dynamically,  inertial,  symmetric  and  convective 
instabilities  are  very  closely  related.  Each  of  these 
instabilities  may  be  thought  of  as  resulting  from  an 
imbalance  of  body  forces  acting  on  a  fluid  parcel.  In  the 
case  of  convective  instability,  the  gravitational  force  acts 
vertically.     In   the   case   of   inertial   instability,   the 
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centrifugal  force  acts  horizontally.  Symmetric  instability 
simply  represents  motion  resulting  from  a  combination  of 
these  forces,  and  thus  the  motion  is  "slantwise."  Following 
Emanuel  (1983a, b),  it  can  be  shown  that  slantwise 
instability  requires  the  slopes  of  the  pseudo  angular 
momentum  (M)  surfaces  to  be  smaller  than  those  of  the 
equivalent  potential  temperature  (  8  e)  surfaces.  This  is 
equivalent  to  stating  that  9e  decreases  along  the  constant 
M  surfaces  upward  or  that  M  decreases  in  the  cross-front 
direction  toward  warm  air  along  a  8e-surface.  From  the 
definition  of  M,  the  latter  implies  that  the  absolute 
vorticity  is  negative  along  0  e-surfaces. 

A  case  study  by  Emanuel  (1983b)  of  a  strongly 
baroclinic  flow  with  observed  linear  bands  of  clouds  and 
precipitation  reveals  that  the  atmosphere  is  initially 
conditionally  unstable  to  upright  convection  and  tends  to 
become  neutral  with  respect  to  slantwise  displacements. 
Such  an  atmosphere  possesses  moist  adiabatic  lapse  rates 
along  M  surfaces,  even  though  it  may  be  substantially  stable 
to  vertical  displacements.  In  this  case,  conditionally 
unstable  baroclinic  flow  undergoes  a  form  of  slantwise  moist 
adiabatic  adjustment  that  renders  the  atmosphere  neutral  to 
slantwise  reversible  displacements.  The  adjustment  uses  up 
the  total  slantwise  convective  available  potential  energy 
(SCAPE)  rather  than  just  the  convective  (gravitational  part 
of)    available    potential    energy    (CAPE)    for    upright 
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convection.  After  such  an  adjustment,  the  SCAPE  of  the 
atmosphere  is  very  close  to  zero  while  the  CAPE  actually 
becomes  negative. 

Emanuel  (1985)  carried  out  a  diagnostic  study  using  a 
semi-geostrophic  model  with  small  potential  vorticity  (i.e., 
small  symmetric  stability)  in  the  region  of  upward  motion 
and  not  so  small  potential  vorticity  in  the  region  of 
downward  motion.  He  demonstrated  that  frontal  forcing 
results  in  a  strong,  concentrated,  sloping  updraft  that  is 
located  ahead  of  the  region  of  maximum  gecstrophic 
f rontogenet ical  forcing.  The  distance  between  the  "front" 
and  the  updraft  will,  in  the  atmosphere,  have  a  typical 
value  on  the  order  of  50-200  km.  This  was  observed  by 
Sanders  and  Bosart  (1985),  who  found  that  slantwise 
convection  can  persist  even  after  a  condition  of  near- 
neutrality  has  been  achieved  under  f rontogenet ical  forcing. 
It  has  also  been  shown  by  Thorpe  and  Emanuel  (1985)  that  in 
an  atmosphere  which  has  small  stability  to  slantwise 
convection,  f rontogenes is  proceeds  at  a  greatly  increased 
rate  (due  to  latent  heat  release)  and  the  horizontal  scale 
of  the  ascent  in  the  warm  air  is  considerably  reduced. 
Cross-frontal  circulations  produced  by  geostrophic  forcing 
are  similar  in  structure  to  slantwise  convection,  despite 
the  (small)  stability  to  the  latter. 
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3 .  Low  Level  Jet  (LLJ) 

In  the  vicinity  of  Japan,  the  Baiu  front  is 
accompanied  by  an  ageostrophic  LLJ  near  700  mb ,  which  is 
oriented  along  the  heavy  rainfall  zone.  The  ageostrophic 
LLJ  is  suggested  to  be  a  result  of  the  downward  transport  of 
horizontal  momentum  caused  by  cumulus  convection  (Akiyama, 
1973a;  Matsumoto,  1972;  Matsumoto  and  Ninomiya,  1969; 
Matsumoto  et  al.,  1971;  Ninomiya,  1971;  Ninomiya  and 
Akiyama,  1974).  In  China,  several  different  relationships 
between  the  low  level  jet  and  heavy  rains  have  been  observed 
in  the  Mei-Yu  season.  In  one  of  these,  a  LLJ  is  coupled  and 
parallel  to  the  upper  tropospheric  jet.  The  formation  and 
strenthening  of  the  LLJ  is  often  associated  with  the  heavy 
rainfall.  Chen  (1982)  suggests  that  coupling  of  upper  and 
lower  level  jets  can  be  explained  by  an  unstable  inertia- 
gravity  wave  that  arises  from  a  thermal  wind  adjustment 
process  that  corrects  an  imbalance  due  to  the  advection  in 
the  entrance  region  of  the  upper  level  jet.  This  coupling 
is  in  contrast  to  that  studied  by  Uccellini  and  Johnson 
(1979)  in  the  Great  Plains  of  the  United  States,  where  the 
LLJ  tends  to  be  perpendicular  to  and  in  the  exit  region  of 
the  upper  level  jet. 

Several  theories  have  also  been  proposed  to  explain 
the  formation  and  diurnal  oscillation  of  the  LLJ  in  the 
Great  Plains  region  of  the  United  States.  Blackadar  (1957) 
suggested   that   oscillation   in   the   boundary-layer   mixing 
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processes  is  responsible  for  the  diurnal  variation  of  the 
LLJ  intensity  which  reaches  maximum  values  by  the  early 
morning.  The  temporal  and  geographic  preference  of  the  LLJ 
have  been  further  explained  by  diurnal  radiation  fluxes  over 
sloped  terrain  (Holton,  1967;  Hsueh,  1970;  Lettau,  1967), 
and  topographical  characteristics  (Paegle  and  Rasch,  1973; 
Paegle,  1978).  Observations  also  showed  that  low  level  jets 
may  also  develop  in  response  to  synoptic  or  subsynoptic 
processes,  particularly  through  a  response  to  leeside 
cyclogenesis  that  is  common  to  the  Great  Plains  (Djuric  and 
Damiani,  1980;  Djuric  and  Ludwig,  1983;  Naistat  and  Young, 
1973;  Newton,  1956;  1967;  Reiter,  1969;  Uccellini,  1980). 
Numerical  simulations  and  a  case  study  by  Uccellini  and 
Johnson  (1979)  suggested  that  the  development  of  the  LLJ  is 
coupled  to  the  upper  tropospheric  jet  streak  by  a  mass 
adjustment  process  within  the  exit  region  of  the  streak. 

C.   OBJECTIVE  OF  THIS  STUDY 

The  objective  of  this  study  is  to  use  a  two-dimensional 
f rontogenesis  model  similar  to  that  of  Williams  et  al. 
(1981)  to  study  the  quasi-steady  state  structure  of  the  Mei- 
Yu  front  and  the  LLJ.  The  model  contains  a  convective- 
adjustment  type  cumulus  heating  parameterization  and  a  bulk 
aerodynamics  formula  for  the  boundary  layer.  To  study  the 
possibility  of  cumulus  momentum  transport  or  the  diurnal 
oscillation  processes,  additional  parameter izat ions  would  be 
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necessary.  Since  considerable  uncertainties  exist  in  both 
the  cumulus  momentum  transport  and  the  diurnal  radiation 
parameter izat ions ,  it  is  desirable  to  carry  out  a  modeling 
study  excluding  these  processes  to  see  if  the  LLJ  may  be 
developed.  One  particular  motivation  for  doing  this  is  the 
tropical  nature  of  the  Mei-Yu  front.  The  intense  convection, 
which  is  correlated  with  the  occurrence  of  strong  low  level 
vorticity,  points  to  a  different  possibility.  The  large 
latent  heat  release  provides  an  energy  source  which  may,  in 
certain  situations,  be  converted  to  the  kinetic  energy  of 
the  LLJ.  This  possibility  is  worth  investigation  especially 
in  light  of  the  recent  symmetric  instability  studies,  which 
show  that  slantwise  convection  may  occur  in  a  frontal  region 
even  if  the  symmetric  instability  condition  is  neutral  or 
slightly  stable. 

Therefore,   the   following   topics   are   pursued   in   this 
study: 

1)  Study  of  the  structure  and  dynamics  of  a  steady 
front,  which  is  investigated  by  varying  the  vertical 
stability,  moisture  content,  Coriolis  parameter  and 
boundary-layer  fluxes  in  the  model;  and 

2)  Study  of  the  relationship  between  the  occurrence  of 
the  LLJ  and  the  convective  activity  in  a  Mei-Yu 
front.  In  particular,  the  possible  role  of  a 
secondary  circulation  equatorward  of  the  front  caused 
by  convection  will  be  explored. 
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In  Chapter  II,  the  basic  formulation  of  the  model  is 
presented.  In  Chapter  III,  the  results  of  the  numerical 
experiments  are  described,  but  the  main  interpretation  of 
the  results  with  respect  to  the  dynamics  of  the  LLJ  is 
discussed  in  Chapter  IV.  Chapter  V  gives  the  summary  and 
conclusions. 


27 


2<n 
22 

20 
18 

16 

MH 

12 

10 

e 

6 

i-\ 

2 


■467A9 


o 

V) 

ea 

QJ 

to 

3 

>- 


ISSi.aiifi 


Q9  ^-^ftOiw=#^wwi 


|05 


06, 


■("■■'i1 — i1'  :  "t    r — i — r — i — r — i — "i •"  i — rr  "i 

Oaie     15    31     15   29    15    31     15    30    15    31     15    30    15    31     15    31     15    30    15    31     15    30    15    31 
Month  1  23  <.  5  676  9  10  11  12 


Fig.   1.1  Mean  daily   rainfall    (y-axis,   mm)    for   the   period   of   1951-1970 

at    Taichung,    Taiwan,    which    is    a    typical    station    having    a 
Mei-Yu  rainfall  maximum  (Chen,  1977). 
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Fig.   1.2  Mean  onset    date   of   the  summer   monsoon   and    the  onset   date 

of  Mei-Yu   season  in  East  Asia   (Tao  and  Chen,   1986). 
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Fig.   1.3  Frequency    of    frontogenesis    at    1°    longitude    by    1°    latitude 

grid  intervals  during  the   Taiwan  Mei-Yu  period  of  May   15  to 
June   15,   1968-1977   (Chen  and  Chi,   1980). 
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Fig.   1.4  Vertical   cross    section    of    (a)    east-west    wind    (u;    solid)    and 

north-south  wind  (v;  dash)  components  in  m  s-1;  (b)  total 
wind  speed  (solid)  in  m  s_1  and  mixing  ratio  (dash)  in 
g  kg-1  along  120°E  at  1200  GMT  11  June  1975.  J  indices 
from  left  to  right  increase  from  south  to  north  with  a  grid 
distance  of  240  km.  The  Mei-Yu  front  at  850mb  is  near  3-1 
(Chen,  1978). 
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Fig.   1.5  GMS-1    satellite    IR    picture    at  0000    GMT    28    May    1981.    a,    b, 

and  c  indicate  organized  mesoscale  convective  system. 
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II.   NUMERICAL  MODEL 


This  study  uses  a  modified  version  of  the  model 
developed  by  Williams  (1974)  and  Williams  et  al.  (1981). 
The  f rontogenesis  in  the  model  is  forced  by  a  nondivergent 
basic  wind  field  that  contains  stretching  deformation  and 
includes  an  Ekman  layer  structure.  The  model  employs 
vertical  and  horizontal  diffusions  of  potential  temperature, 
moisture  and  velocity.  The  diffusion  coefficients  vary  only 
with  elevation.  Pressure  coordinates  are  used  in  the 
vertical  and  a  tropopause  is  simulated.  This  should  lead  to 
a  more  accurate  simulation  of  the  upper  regions  where  latent 
heating  will  be  important.  The  time-dependent  quantities  in 
the  model  are  functions  of  y  (cross-front)  and  p  and  the 
model  is  integrated  until  a  quasi-steady  state  is  achieved. 
A  planetary  boundary  layer  is  added  with  aerodynamic  bulk 
formulae  for  the  surface  fluxes.  A  moist  convective 
adjustment  process  is  included  to  represent  cumulus 
convection.  A  staggered  grid  is  used  in  the  horizontal 
because  it  gives  a  better  treatment  of  the  solutions  that 
involveheating. 

In  Section  A,  the  basic  equations  in  pressure 
coordinates  are  developed  including  the  appropriate  moisture 
terms.  The  convective  adjustment  and  latent  heating  are 
discussed  in  Section  B.  In  Section  C,  the  boundary  layer 
treatment   is   described.      The   initial   conditions   and 

33 


tropopause  simulation  are  in  Section  D.  The  details  of  the 
specification  of  model  grids,  grid  arrangements  and  finite 
difference  scheme  are  in  Appendix  A,  B  and  C  respectively. 

A.   BASIC  EQUATIONS 

The    hydrostatic    primitive    equations    in    pressure 

coordinates  may  be  written 


av -        (2.1.1) 

~  +  v  ^vvj   +  -y-  udv;    t  vq>   t    nun   -  ^mv     v    ■     -5p     L  v-^«"5p- 


8V  +V-(VV)    +  -|r:  M)    +  v<»   +  fkxV  =  Amy2V  +^7    [WP§£1    • 


||  +  v  .(ev)  +  -^(Qu)    =  Aev2e  +  -^  [k(p)  |L]  +  oA  +  Qm  +  hc  , 

|f +  V.(qV)   +-|p   G»q)   =Aqv2q+^[K(p)^]  +  MA-Mc    ,  (2.1.3) 


V  .V  +  5co/op   =  0, 


(2. 1.4) 


54>/op  =  -r(p)6,  (2.1.5) 

where: 

V  =   horizontal    component   of  the  velocity, 

to  =  dp/dt, 

p  =  pressure, 

$  =  gZ  =  geopotential , 

e  =  T(Po/p)<    =  potential    temperature, 

T  =  temperature, 

R  =  gas   constant, 

Cn  =  specific  heat  at  a  constant  pressure, 

k  =  R/Cp, 
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f  =  2Q    sin  cj;    =  Coriolis   parameter, 

Q   =  rotation   rate  of  earth, 

<\>   =  latitude, 

q   =  specific   humidity, 

Am  =  horizontal    diffusion  coefficient   of  momentum, 
Aq    =  horizontal    diffusion  coefficient   of  heat, 
A«  =  horizontal    diffusion   coefficient   of  humidity, 
K(p)   =  vertical    diffusion  coefficient, 
Q/\  =  heating  due  to  dry  convective  adjustment, 
Ofj\  =  heating  due  to  moist  convective  adjustment, 
H£  =  large  scale  latent  heating, 

Mfl  =  specific  humidity   change  due  to  moist   convective  adjustment, 
Mq  =  specific  humidity  change  due  to  large-scale  condensation, 

r(p)  -Rp-i (_£_,*.  cp-^f-jj-*. 
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The  upper  and  lower  boundary  conditions  are  given  by 


w(x,  y,  pt,  t)  =oj(x,  y,  p0,  t)  =  0, 
V(x,  y,  p0,  t)  =  0, 

§-p(x,  y,  Pt,  t)  =  0, 

^p(x,  y,  pt,  t)  =  g-jj-(x,  y,  p0,  t)  =  0. 


(2. 1.6) 


Here  pt  and  po  represent  the  pressure  at  the  top  and  bottom 
of  the  atmosphere,  respectively,  and  these  quantities  are 
approximated  as  constants.  The  conditions  (2.1.6)  include 
the  kinematic  boundary  condition  at  the  top  and  bottom  of 
the  atmosphere,  the  zero  slip  condition  at  the  surface,  the 
zero  stress  condition  at  the  top  of  the  atmosphere,  and  the 
zero  heat  flux  conditions  at  the  top  and  bottom  of  the 
atmosphere,  respectively. 

The  approximate  steady-state  solutions  to  (2.1.1)- 
(2.1.5)  that  satisfy  the  boundary  conditions  (2.1.6)  are  as 
f ol lows : 

Ti  =  D[x(l  -  e<   cos  5)  +  ye<    sin  £l 


v 

U) 

0 

q 


=  -n[y(l  -  e<   cos  I)   +  xe<    sin  &] 

=  0, 

■S  s  -D2(x2  +  y2)/?.   -  fDxy  +F(p), 

=  R(p). 

■  H(p), 


(2.1.7) 
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where  £ =  [ f /2K( po ) ] l / 2 ( po  -  p).  The  deformation  D,  and  the 
Coriolis  parameter  f  are  constant,  and  the  functions  F(p), 
G(p)  and  H(p)  are  arbitrary  except  that  F  and  G  must  satisfy 
the  hydrostatic  equation  (2.1.5).  These  solutions  are  only 
approximate  within  the  Ekman  layer  because  some  small 
advection  terms  have  been  neglected. 

The  dependent  variables  are  subdivided  as  follows: 


U  =  u(x,  y,  p)  +  u(y,  P,  t), 

V  =~v(x,  y,  p)  +  v(y,  P,  t), 

(D  =  w(y,  p,  t),  9  =  9  (y,  p,  t), 

<t>  =$(x,  y)  +-rc(y,  p,  t),  q  =  q(y,  p,  t), 


(2. 1.8) 


where  all  of  the  departures  from  the  steady-state  solutions 

are  assumed  to  be   independent  of   x.      If  we  substitute 

the   relations   (2.1.8)   into   (2.1.1)  -   (2.1.5)    and   use 
dv 


5_U 
5x 


=  0(1  -  e<   cos  Fj 


5u    5  v 


=  De~£  sin  Z, ,  u,  v  satisfy 


oy   "x*   -   —  -»#•    5y    5x 
the  steady-state  equations.   The  equations  of  the  departure 

field  from  the  steady-state  solution  become 


5_u 
ot 


|_(uv)+  4?(w)+u|i  +  -,|i  +  v^  +  uau 


op 


=   fv   +  A 


52u 
n  5^- 


5      (v  aux 

Tp   (K  oT}    * 


(2.1.9) 
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St  +  "ay  (vv)  +  Tp  (wv)  +  Ty  (vv)  +  u  5~  +  w  aT 


.      bn      .    a2v         a     fv  avv 
■  "fu  "  ay  +  A*  W       "^P  {K  a?} 


(2,1.10) 


ae       a     /oN       a     ,AN   L  ~  ae      .    o2e       a     ft/  aeN 
at  +  "ay  (ve }  +  Tp  M }  +  v  57  =  **  &F     Tp  (K  aT} 


+  0A  +  0M  +  HC, 


(2.1.11) 


+  MA  -  MC    , 


(2.1. 12) 


—    +    —        n 

ay         aP 


(2. 1. 13) 


a^         /  \« 
aT  =  "r(p)e 


(2. 1. 14) 


Using  the  relationship  of  (2.1.7),  the  (2.1.9)  becomes 


a  11  ,    a     /    »       5     ,  .  - 

•  (uv)   +  -^—  (uu)   +  0(1   -  e<    cos  l)u 


at      ay 


-  n[y(l  -e<   cos  I)   -  xe<   sin  ll~+  D{e<   sin  Ov  (2.1.15) 


9y 


+  wn  Tp"  Wl   "  e"^   cos  ^)  +  ye"?    sin  5  1 


=  fv  +  Am  L4"  +  4-  (K  P-)    . 


Am 


ar       ap  v^  ap 
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It  is  to  be  noticed  that  the  assumption  that  u  is 
independent  of  x  for  all  time  is  violated  in  (2.1.15)  in  the 
two  terms  where  x  appears.  This  indicates  that  the  other 
quantities  which  were  assumed  independent  of  x  are  actually 
x-dependent,  which  would  bring  other  terms  into  the 
equations.  In  this  study,  these  added  terms  are  neglected 
and  (2.1.15)  and  the  other  equations  are  applied  at  x=0.  If 
the  initial  u,  v,  &  ,  9  and  q  fields  are  independent  of  x, 
the  error  from  this  approximation  will  grow  slowly  and  it 
will  be  confined  to  the  boundary  layer  since  the  x-dependent 
terms  in  (2.1.15)  are  zero  outside  the  boundary  layer. 
Since  frontogenesis  occurs  very  rapidly,  it  is  expected  that 
the  development  of  x-variations  in  the  dependent  variables 
would  only  have  a  small  effect  on  the  resulting  quasi-steady 
front.  In  any  case,  these  effects  could  not  be  observed  in 
the  atmosphere,  because  atmospheric  fronts  always  have  some 
variation  in  the  basic  fields  along  the  front  which  would  be 
much  more  important. 

When  (2.1.15)  is  evaluated  at  x  =  0 ,  it  becomes 


+  \){e<    sin  I  )v  -  wVn  —  (e<    sin  l)   =  fv  +  A 


dp 


Am 


dT"  Tp"  {K  a~p' 


(2. 1. 16) 


39 


where   Vq  =  -fly   and  T    =  V^(l  -  e~^  COS  £).       To   remove   the 
vertical  mean  in  (2.1.10),   the  vertical  average  is  defined 


as 


<(    )>  = 


Po  -  Pt 


/   (   )dp. 
Pt 


By  vertically  averaging  (2.1.13)  and  using  the  boundary 
conditions  (2.1.6),  we  obtain 

<V>  =  n  . 


ay 


This  equation  states  that  the  total  disturbance  y  mass  flux 
is  independent  of  y.  It  can  be  expected  that  at  a  large 
distance  from  the  frontal  zone  the  disturbance  mass  flux 
will  vanish.   Thus,  we  get 


<V>  =  0  . 


(2.1.17) 


If  we  take  the  vertical  average  of  (2.1.10)  and  use  (2.1.6) 
and  (2.1.7),  we  obtain 


_L_  <vv>  +  4-  <vv>  +  <u  U>  +  4o  ¥r> 
dy  5y 


ox 


dp' 


B*>   .   f<u>   +   [JL|^] 
dy  Po  5P 


Po 


(2.1. 18) 
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When    this    result    is    subtracted    from    (2.1.10),    we    have 


|v  +  _|_  (vv  .  <vv>)  +  4_  (MV)  +  -§_  fvv  -  <vv» 


ov       y    5v,  x      ov       ,    5v. 

+  ,J  5T  "  <lJ  Sx*  +  w  5T  "  ^  a?" 


=     -    f(u     -     <V>)     -   -|-    (lC     -     <TC>) 


5^v 


o      /  w  a  v  n        ,   K    '5  v . 


+   V  g^f  +       5p    (K  dp^    "    (    p0  op 


Po 


Using      (2.1.7)    and    evaluating    at    x    =    0,    we    obtain 


U  +  T7  (vv  "  <vv>)  +  4p  (uv)  +  T?  (rv  "  <rv>) 


+  n[ue"£   sin  £,   -  <ue"^   sin  £>  -  Vn[< 


00 


op 


{&<   cos  g)   -  <fc>~-     (e<   cos  £)>] 


|y    («     -    <1X>)     -    f(U     -    <V>)     +    Am|^ 


+    aP  (K5P}      (  PoaP) 


(2.1. 19) 


Po 


Eq.(2.1.11)    evaluated    at    x    =    0    may    be    written, 


d9 
5 


t        ay        apv;        ay 


(2.1.20) 


a2e 


ae 


Ae  ^  +  -^rK(P>oT]  +  °a  +  Om  +  hc  . 
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Similarly,     (2.1.12)    becomes 


St+-S7(v"i+Tp  W& 


=  Aq|i^  +  _|^  :k(p)  |^n  +  ma  -  Mc 


(2.1.21) 


Integrating'   the   hydrostatic   equation   (2.1.14)   with 
respect  to  p  and  removing  the  vertical  mean,  we  obtain 


%   -  <rc>  =  J   r(p)9dp  -  <J   r(p)0dp  >  . 
Pt  Pt 


(2. 1.22) 


With  (2.1.13),  (2.1.16)  and  ( 2 . 1 . 19) - ( 2 . 1 . 22 )  form  a 
complete  set  of  equations  which  can  be  solved  by  a  marching 
process . 


R.   CONVECTIVE  ADJUSTMENT  PROCESSES  AND  LATENT  HEATING 

This  section  describes  the  dry  and  moist  convective 
adjustment  processes  and  the  latent  heat  release  that  are 
used  in  the  model.  The  dry  convective  adjustment  process  is 
activated  in  any  layer  where  89  /9p  >  0.  In  such  a  layer, 
6  is  set  equal  to  the  average  value  within  the  layer.  Then 
the  layers  above  and  below  are  tested  for  instability  and 
modified  if  necessary.  This  process  is  represented  by  Qa  in 
(2.  1.20)  . 

When   the   relative   humidity   exceeds   100%,   that   is 
q  >  qs  [qs(T,p)  is  the  saturation  mixing  ratio],  latent  heat 
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of  condensation  will  be  released.    To  adjust  the  relative 
humidity  to  100%,  T  and  q  must  be  changed  to  satisfy 


q  +  5q  =  qs(T  +  6T,  p) , 


Cp6T  =  -L5q, 


(2.2.1) 


(2.2.2) 


where  L  is  the  latent  heat  of  condensation,  and  6  <1  and  <$  T 
are  the  changes  in  q  and  T  during  the  process.  A  good 
approximation  to  the  function  qs(T,p)  is  given  by 


Qs 


H.38     r  L  ,        1 
-T-   exp[~R^  (  273°K 


)], 


(2.2.3) 


where  Rv  is  the  gas  constant  for  water  vapor. 

The  equations  (2.2.1)  and  (2.2.2)  require  that  the  final 
state  of  the  air  be  exactly  saturated  and  that  excess 
moisture  be  condensed  isobarically ,  which  releases  the 
latent  heat  to  the  air. 

Equations  (2.2.1)  and  (2.2.2)  can  be  expressed  in 
iterative  form  as 


6q  =  qs  (1  -  r)  +  (-|^S-)  6T,   q  =  rqs 


(2.2.4) 


and 


cp6T  =  -L6q 


(2.2.5) 


Eliminating   <5  T  from  (2.2.3)  and  (2.2.4)  gives 
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6q  = 


cp  l  5T   ;p 


(2.2.6) 


At  constant  pressure,  the  C lausius-Clapeyron  equation  can  be 
expressed  as 


aT  'p 


"R7  ~V 


(2.2.7) 


Substituting  (2.2.6)  into  (2.2.7)  gives 


6q  = 


9sB  -  r) 


l  + 


L2 


3s 


CpRyT 


(2.2.8) 


<5  q  may  now  be  calculated  from  (2.2.8)  as  can  the 
corresponding  value  of  6T  from  (2.2.4).  The  first  estimates 
of  the  adjusted  temperature  and  mixing  ratio  are, 


T'  =  T  +  6T    and    q'  =  q  +  6q  . 


(2.2.9) 


These  values  are  in  turn  used  to  obtain  improved 
approximations.  Iterations  are  continued  until  the  values 
of  T'  and  q'  at  each  level  match  with  sufficient  accuracy  a 
value  of  100%  relative  humidity.  The  changes  are  then 
applied  to  the  potential  temperature  and  mixing  ratio  fields 
of  the  model,  which  accounts  for  the  term  He  of  (2.1.20) 
and  Mc  of  (2.1.21). 
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The  convective  latent  heating  is  parameterized  by  a 
convective  adjustment  scheme.  The  scheme  yields  similar 
results  in  describing'  the  large-scale  effect  of  latent 
heating  as  an  explicit  scheme  which  included  prediction 
equations  for  water  vapor,  cloud  water  and  rain  water  (Hsie 
and  Anthes,  1984) . 

In  areas  where  convection  is  occurring,  it  is  necessary 
to  determine  whether  the  atmosphere  is  convectively  unstable 
to  moist  processes.  A  layer  will  be  conditionally  unstable 
if  30e / 9p  >  0,  where  0e  is  the  equivalent  potential 
temperature,  which  is  given  by 


ee  =  e  exp(-^_) 


(2.2. 10) 


As  with  the  dry  convective  adjustment,  6e  is  set  equal  to 
the  average  value  in  the  unstable  layer.  This  gives  a  value 
for  9e  at  each  level  and  6T  and  6q  can  then  be 
determined  from  the  following  calculations. 

The  heat  release,  Qm  ,  as  a  result  of  moist  convective 
adjustment  is  included  through  knowledge  of  the  change  of 
equivalent  potential  temperature  during  the  adjustment. 
Taking  the  natural  logarithm  of  (2.2.10)  gives 


in  9e  =  In   9  +  ■  Lc^ 


(2.2. 11) 
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A  constant  pressure  surface  may  be  assumed  since  the  model 
uses  pressure  coordinates.  Differentiating  both  sides  of 
(2.2.11)  and  recalling  the  definition  of  potential 
temperature  gives 


66, 


6T 


6q<;L 
cpT 


Lq«£T 

T2- 


(2.2. 12) 


Solving  (2.2.6)  for  dqs ,  substituting  into  equation  (2.2.11) 
and  solving  for  dT  leads  to 


6T  = 


69 


eef[(i 


3s 


RyCpT 


L2     Lq, 


cpT 


(2.2.13) 


where    <$0e   is   the   change   in  6e   needed   to   adjust   the 

gradient  30e/3p  to  be  equal  to  zero. 

Since   the   atmosphere   remains   saturated  throughout   the 

process   of   moist   convective   adjustment,   the   Clausius- 

Clapeyron  equation  is  used  again  to  obtain  the  corresponding 

6q   or   actual   5 qs .     Changes   are   again   applied   to   the 

potential  temperature  and  mixing  ratio  fields  of  the  model 

which  accounts  for  the  terms  Qm  of  equation  (2.1.20)  and  Mc 

of  equation   (2.1.21).     The  convective   adjustment   with   a 

saturated  atmosphere  is  a  rather  crude  parameterization  of 

the  cumulus  effects.    However,   the  relative  humidity  can 

indeed  become  very  high  in  intense  Mei-Yu  cases   (Chen  and 

Tsay,  1978).   Furthermore,  the  strong  control  by  large-scale 
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frontogenesis  forcing  makes  the  results   less  sensitive  to 
the  details  of  the  cumulus  parameterization. 

C.   BOUNDARY  LAYER  PARAMETERIZATION 

The  surface  fluxes  are  parameterized  by  the  bulk- 
aerodynamic  formula  and  the  surface  fluxes  of  heat  and 
momentum  are  assumed  to  influence  directly  only  the  lowest 
model  layer.  The  magnitude  of  the  surface  stress  is 
proportional  to  the  square  of  the  wind  speed  and  the  surface 
sensible  heat  flux  is  proportional  to  the  product  of  the 
wind  speed  and  the  vertical  gradient  of  temperature  between 
the  ground  and  the  middle  of  the  lowest  model  layer.  This 
formulation  is  best  suited  to  models  with  coarse  vertical 
resolution,  since  the  depth  of  the  lowest  layer  should 
correspond  to  a  typical  PBL  depth. 

In  this  study,  the  vertical  sub-grid  scale  mixing  terms 
in  Eqs. (2.2.16)  are  applied  to  the  lowest  layer  of  the 
model.   The  term  may  be  written  as 


hiKMU]  =VF"(p)] 


(2.3.1) 


At  the  top  of  the  boundary  layer,  we  assume 


Fu(p)  =  Fu(z)  =  K(z) 


5_U 

oz  * 


(2.3.2) 


K(z)  is  defined  as 


K(x)  =  CD|u| 
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here  Cd  is  the  drag"  coefficient  and  |  u  |  =  +  /  u2  +  v2 
Applying  (2.3.2)  to  the  lowest  two  layers,  po  and  pi,  of  the 
model ,  we  get 

Fu  =  "  P9  Cf)|u|(ui  "  U0)/(P1  "  PO) 


-  P9  CD I  u  I  u  1  /A  p 


since  uo  =  0  from  the  lower  boundary  condition  and  Ap  =  pj  -  PQ  . 
Therefore,  the  u  momentum  flux  at  the  surface  becomes 


Fu  .  .  in-3  (    PQCn  |u|  }  U1  t 


(2.3.3) 


where  10~3   is  a  unit  conversion  factor  which  converts  the 
unit  of  pressure  to  cb .   Similarly,  the  v  momentum  flux  is 


Fv_=  -  10-3  (pgcp  H  )vi  ^ 


(2.3.4) 


In  this  study,  the  PBL  depth,  h,  is  estimated  as 


h   = 


po  -   Pl 
pg 


4^S  x  *  ■ « -■ 


The  sensible  heat  and  moisture  fluxes  are 


Fe=.in-s  (Ps;-nH  )  (6l-e0) 


Fq  = 


-  m-3  ("ft1"1  )  <«,  -  %'• 


(2.3.5) 
(2.3.6) 
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where  9o  and  qo  are  the  specified  surface  values 
The  latent  heat  flux  is  defined  as 


FL  ■ 


L  F, 


(2.3.7) 


where  L  =  2.5  x  106  Joule/kg,  cP  =  1003  Joule/ (° K/kg) 


D.   INITIAL  CONDITIONS 

The  initial  potential  temperature  is  given  by 


e(y»  P>  n)  =  G  i  (p)  -  a(-~—)  arctan  (sinh  ay)  , 


(2.4. 1) 


where  8  i  (p)  is  a  prescribed  vertical  potential  temperature 
profile  and  a  is  one-half  the  total  horizontal  potential 
temperature  variation.    a  is  defined  as  and  is  evaluated 

at  some  middle  level  p  =  pm . 


*/h/X  (Xf  (-Mi) 

Pm   Pm      ^P 


Pm 


The  thermal  wind  equation,  which  is  derived  using  the 
hydrostatic  equation  and  the  definition  of  geostrophic  wind, 
may  be  written  as 


&U.  _  1  R.  /_p_xK  59_ 
5P    fp1  P0;  Sy  ' 


(2.4.2) 


Integrating  (2.4.2)  from  p  =  pm  to  p  -    p 


pis~dp"78~  p!  p(^  dp 


and  assuming  u  =  0  at  p  -  pin,  we  obtain 
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u(,.p.o).  %  |i[(j-,<  .(^n 


Po  Po 


Using    (2.4.1), 
59_ 

ay 


2a    5 


%     dy 
2aa 


farctan    (sinh  ay)1 


% 


sech  (ay)  . 


(2.4.3) 


(2.4.4) 


Substituting  (2.4.4)  into  (2.4.3),  the  initial  u  component 
of  the  velocity  may  be  written, 


2axC 


u(y,  p,  0)  =  -  ^P  [ML-)K  -  (£mf  1  sech  (ay) 


Po 


Po 


(2.4.5) 


The  v  component  of  the  velocity  is  set  initially  to  zero, 


i.e.  , 


v(y,  p,  o)  =  o  , 


so  that ,    to  (y,  p,  0)    =  0  . 

The  initial  field  of  mixing  ratio  is  calculated  by 
determining  the  saturation  vapor  pressure  es  at  each  grid 
point  from  an  integrated  form  of  the  Clausius-Clapeyron 
equation.   According  to  Hess  (1959), 


es  =  0.611  exp  (-p—  (— —  j  ))     , 


(2.4.6) 


where   To   is   a   constant   of   integration.     Knowing   the 
saturation  vapor  pressure,   the  saturation  mixing  ratio   is 
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obtained  from  the  approximate  relation  (Haltiner  and  Martin, 
1957) . 


qs  = 


0.622  e< 


(2.4.7) 


Knowing  the  potential  temperature  at  t  =  0  from  (2.4.1),  it 
is  possible  to  calculate  qs  at  each  point  using  equations 
(2.4.6)  and  (2.4.7).  Defining  q  to  be  a  given  fraction  of 
qs  allows  the  testing  of  atmospheres  with  different  initial 
relative  humidities. 

The  initial  conditions  described  in  previous  sections 
will  be  constrained  by  a  simulated  tropopause  that  is 
centered  at  about  300  mb  level.  The  simulation  is 
accomplished  by  using  a  reference  temperature,  Tref ,  at 
about  the  midpoint  of  the  horizontal  scale  and  calculating 
temperatures  such  that  a  tropopause  with  the  expected 
variation  in  height  is  determined  and  a  zero  vertical 
temperature  gradient  is  achieved  in  the  lower  stratosphere. 
A  reference  temperature,  Tref,  is  computed  using  (2.4.1) 


Tref  =  215  +  a  (-J-)  arctan  (sinh  ay)  (  Pn  -  Prpf f    .    (2.4.8) 

n  Po 


When  the  temperature  in  each  column  reaches  this  reference 
temperature,  it  is  held  constant  at  that  value,  as  expected 
in  the  lower  stratosphere.  The  temperature  is  then 
converted  to  potential  temperature. 
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III.   NUMERICAL  EXPERIMENTS  AND  RESULTS 

In  this  chapter  two  groups  of  numerical  experiments  are 
carried  out  to  simulate  the  midlat i tude  and  subtropical 
situations,  respectively.  For  each  group,  the  Coriolis 
parameter  and  initial  temperature  sounding  are  fixed,  but 
the  amount  of  the  initial  relative  humidity  and  surface 
fluxes  of  momentum,  mixing  ratio  and  sensible  heat  are 
varied.  The  purpose  here  is  to  describe  the  structure 
differences  due  to  the  difference  in  latitude,  relative 
humidity  and  surface  fluxes.  The  discussion  of  the  results, 
which  focuses  on  the  development  of  the  low  level  jet,  is 
given  in  Chapter  IV. 

We  will  first  give  the  parameter  ranges  used  in  the 
experiments  (in  Section  A),  and  then  describe  the  numerical 
results  of  the  midlatitude  (Section  B)  and  subtropical 
(Section  C)  cases,  respectively.  As  mentioned  in  the 
Introduction,  the  main  purpose  for  all  numerical  experiments 
is  to  obtain  a  quasi-steady  solution.  However,  it  is  also  of 
interest  to  follow  the  time  evolution  of  the  f rontogenesis 
process  to  see  how  the  quasi-steady  states  are  developed. 
Therefore,  in  each  group  the  evolution  for  a  "typical  case" 
is  shown.  It  turns  out  that  these  evolutions  provide  useful 
insights  about  the  different  nature  of  the  effect  of  cumulus 
convection  in  the  midlatitude  and  subtropical  cases. 
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A.   EXPERIMENT  CASES  AND  PARAMETERS 

The  following  values  are  used  for  the  constants  in  this 
study : 


9 

=  9.81  m  s-2  . 

cp 

=  1003  joule/kg/°K 

R 

=  287  joules/kg/°K 

K 

=  R/Cp 

Rv 

=  461  joules/kg/°K 

L 

=  2.5  x  106  joules/kg 

57 

=  5°K/km 

(59j 

[  ap 

)Pm 

=  -  0.53448°K/cb 

Po 

=  100  cb 

Pm 

=  55  cb 

% 

=  Aq  =  %  =  3  x  104  m2  /s 

Km 

=  Kq  =  Kq  =  0.0007  cb2/s 

At 

=  270  sec 

Ay 

=  60  km 

Y 

=  2760  km 
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Each  experiment   is   designated  by   a  three-letter   case 
name  as  follows : 

LRF 
where  L  defines  the  latitude  environment. 

1.  L  =  M  refers  to  a  midlatitude  environment,  with  the 
following  parameter  settings: 

a.  f  =  1.06  x  10-4  s"1  (corresponds  to  46°N); 

b.  initial  vertical  potential  temperature 
distribution  given  by  a  typical  midlatitude 
sounding  (Fig.  3.1); 

c.  initial  meridional  temperature  variation  is 
given  by  a  =  8°K;  and 

d.  large-scale  deformation  fields  given  by  D  = 
1  x  lO"5  s"1  . 

2.  L  =  S  refers  to  a  subtropical  environment,  with  the 
following  parameter  settings: 

a.  f  =  0.558  x  10~4  s"1  (corresponds  to  22.5°N); 

b.  initial  vertical  potential  temperature 
distribution  given  by  a  typical  subtropical 
sounding  (Fig.  3.1); 

c.  initial  meridional  temperature  variation  is 
given  by  a  -  4°K;  and 

d.  large-scale  deformation  fields  given  by  D  = 
0.5  x  10-5  s-1  . 
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3.  R  defines  the  relative  humidity: 

a.  R  =  6  means  R.H. =60%;  and 

b.  R  =  9  means  R.H. =90%. 

The  observed  low-level  relative  humidity  in  the  vicinity 
of  a  Mei-Yu  front  varies  between  60%  -  90%  (Chen  and  Tsay, 
1978).  When  the  Mei-Yu  front  is  north  of  the  Yantze  River  (a 
midlatitude  situation)  the  minimum  R.H.  is  about  60%.  When 
it  is  located  over  Taiwan  and  southeastern  coast  of  the 
Chinese  mainland  the  maximum  R.H.  is  near  90%. 

4.  F  defines  the  surface  fluxes: 

a.  F  =  0  means  there  is  no  surface  fluxes  of  heat, 
moisture,  and  momentum;  and 

b.  F  =  1  means  surface  fluxes  of  heat,  moisture  and 
momentum  are  included. 

Typically,  the  fluxes  are  large  when  a  subtropical 
portion  of  the  front  is  over  a  warm  sea  surface.  When  the 
front  is  moved  over  land  in  the  midlatitude  cases,  the 
fluxes  are  normally  weaker  with  the  exception  of  momentum 
fluxes.  For  self-consistency,  the  fluxes  of  all  quantities 
will  either  be  included  or  excluded  in  the  numerical 
experiments.  For  all  cases,  the  initial  zonal  wind  is 
calculated  from  the  potential  temperature  distribution  using 
the  thermal  wind,  and  the  initial  v  component  of  the 
velocity  is  set  to  zero. 
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B.   MIDLATITUDE  SIMULATIONS 

1 .   Typical  Midlat i tude  Case  (M60) 

The  midlatitude  case  of  M60  with  an  initial  R.H.=60% 
and  no  surface  fluxes  describes  the  development  of  the  front 
in  detail.  The  potential  temperature  and  mixing  ratio  at 
t=0  are  shown  in  Fig.  3.2a,  and  the  zonal  velocity  at  t=0  is 
shown  in  Fig.  3.2b. 

Fig.  3.3  depicts  the  development  of  the  potential 
temperature  and  mixing  ratio  fields.  The  horizontal 
gradient  of  both  fields  increases  steadily  through  the  first 
two  to  three  days  as  a  frontal  structure  develops  near  42°- 
43° .  The  maximum  surface  temperature  gradient  nearly 
triples  (from  the  2°K/degree  latitude  to  about  6°K/degree 
latitude),  while  the  maximum  surface  mixing  ratio  gradient 
increases  nearly  ten-fold  (from  about  0.5  g/kg/degree 
latitude  to  about  5  g/kg/degree  latitude).  This  midlatitude 
mixing  ratio  gradient  is  about  twice  as  strong  as  that 
observed  in  a  subtropical  Mei-Yu  front  by  Chen  (1978)  as 
shown  in  Fig.  1.4b,  and  the  average  mixing  ratio  is  lower 
(8  g/kg  versus  11  g/kg)  .  The  axis  of  maximum  mixing  ratio 
exhibits  a  conspicuous  poleward  vertical  tilt  from  41°  at 
the  surface  to  48°  at  8  km.  This  poleward  vertical  tilt 
with  a  slope  of  approximately  1.1:100  agrees  well  with  the 
typical  midlatitude  fronts  (Byers,  1974).  The  tilt  is 
slightly  stronger  than  that  observed  in  the  subtropical  Mei- 
Yu   front   (slope   1.3:100)   as   shown   in   Fig.   1.4b.     The 
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development  of  the  frontal  thermal  structure,  including  the 
concentration  of  the  isentropes  in  the  lowest  2  km  and  the 
poleward  tilt  of  the  frontal  zone  from  the  surface  to  the 
upper  troposphere,  is  clearly  indicated  in  Fig.  3.4.  In 
addition,  a  secondary  maximum  occurs  in  the  upper 
troposphere  (6  -  8  km  or  600  -  400  mb )  near  50°,  and  a 
reversed  tilt  and  reversed  gradient  is  present  in  the  lower 
stratosphere.  This  stratospheric  structure  is  a  result  of 
the  imposed  condition  of  a  tropopause  and  is  consistent  with 
previous  results  obtained  by  Williams  et  al.  (1981). 

The  vertical  velocity  field  (Fig.  3.5)  quickly 
develops  and  becomes  organized,  with  a  concentrated  upward 
motion  of  >■  -4.0  x  10~4  cb/sec  occurring  at  day  2  in  the 
lower  troposphere  near  43°.  A  poleward  vertical  tilt  of  the 
area  of  upward  motion  develops  and  becomes  more  obvious  with 
time.  After  day  4,  a  quasi-steady  state  clearly  delineates 
both  the  upright  maximum  rising  motion  centered  at  43°  to 
44°  in  the  lower-middle  troposphere  and  the  vertical  tilt, 
which  is  well  represented  by  the  (a  -0  line  poleward  of  the 
maximum  rising  motion. 

A  very  small  area  with  a  secondary  maximum  in 
vertical  velocity  appears  at  day  2  near  the  surface  ahead  of 
the  front  (41°).  This  represents  the  effect  of  the  forcing 
by  the  Ekman  layer  structure  of  the  mean  flow,  as  was 
discussed  by  Williams  et  al.  (1981).  Because  this  area  is 
very  close  to  the  surface,   it  does  not  cause  condensation 
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and   latent   heat   release   even   in   the   highly  moist   cases 

discussed  below. 

The  initial  zonal  wind  (Fig.  3.2b)  is  in  geostrophic 

balance   with   the   temperature   field,    with   an   upper 

tropospheric  westerly  -maximum  of  ~55  m/s   and  a   surface 

easterly   maximum   of   -25  m/s.     The   pattern   is   almost 

symmetric  with  respect  to  46°  latitude.   As  the  integration 

proceeds  (Fig.  3.6),  asymmetries  develop  and  a  strong  shear 

begins  to  appear  at  day  1.0  below  6  km.    The  shear  quickly 

reaches  a  steady  state  with  a  frontal  structure  which  tilts 

from  42°  at  the  surface  towards  45°  at  4  km.   This  is  in  the 

same  vicinity  as  the  go  =0  line  and  the  frontal  zone  that  is 

defined   by   the   potential   temperature   and   mixing   ratio 

gradients.    At  the  latitude  of  the  surface  front,  a  local 

westerly  maximum   develops   near   z  =  4  km.    By   day  2   (Fig. 

3.6b),  this  narrow  westerly  regime  penetrates  downward  and 

appears   to   be   an   intrusion   into   the   lower   tropospheric 

easterlies.    It  remains  in  the  3-4  km  layer  through  day  4. 

Although  the  magnitude  of  westerly  wind  in  this  narrow  zone 

is  small,   the  poleward  meridional  shear  is  large,    which 

results  in  a  vorticity  in  this  region  that  is  comparable  to 

the  values  near  the  upper  tropospheric  westerly  jet.    The 

development  of  this  narrow  westerly  regime  is  a  result  of 

the  intensifying  meridional  temperature  gradient  which  is 

forced   by    the    f  rontogenesis    process    in    the    lower 

troposphere.   Development  of   this  zonal   velocity  field   is 
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also  dramatically  illustrated  by  the  evolution  of  the 
relative  vorticity  (Fig.  3.7).  The  initial  flow  pattern 
develops  rapidly  into  a  quasi-steady  state  frontal  structure 
after  about  two  days.  Maximum  cyclonic  vorticity  in  the 
lower  troposphere  is  concentrated  on  the  warm  side  of  the 
front.  To  illustrate  the  vertical  structure  of  the 
horizontal  wind  shear  with  that  of  the  horizontal 
temperature  gradient,  3  6  /  9y  =  -2°K/100  km  for  day  4  is 
plotted  in  the  vorticity  figure.  A  coherent  vertical  tilt 
of  the  frontal  structure  is  obvious. 
2.   Other  Midlatitude  Cases 

The  effects  of  varying  the  initial  relative  humidity 
(M90),  and  the  inclusion  of  surface  momentum,  heat  and 
moisture  fluxes  (M61  and  M91)  are  illustrated  by  the 
resulting  structure  at  day  4  (Figs.  3.8-3.12).  In  addition, 
the  precipitation  and  heating  rates  are  also  shown  in  Figs. 
3.13  and  3.14,  respectively. 

a.   Effect  of  Moisture 

(1)  No  Surface  Flux  Cases  (M60  vs.  M90) 

It  is  clear  that  the  initial  R.H.=90%  case 

(Fig.  3.8b)  has  higher  mixing  ratios  in  general  and  a  warmer 

upper   troposphere   between   40°-50°   (best   represented   by 

comparing  the   9=350°K  line).    The  higher  temperature  is  a 

result  of  a  larger  latent  heat  release  in  the  frontal  region 

due  to  the  increased  moisture  content,   as  can  be  seen  by 

comparing    the    precipitation    rates    (Fig.    3.13).    The 
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meridional   potential   temperature   gradient   in   the   surface 

front  for  the  R.H".=90%  case  is  also  stronger  (Figs.   3.9a 

versus  3.9b).    The  frontal  zones  as  defined  by  the  region 

36/3y   <  -0.02°K   km-1   are   quite   similiar   in   width   and 

vertical  tilt.   Above  the  surface  front  position  (~42°),  the 

mixing  ratio  in  the  R.H.=90%  case  (Fig.  3.8b)  also  shows  a 

significant  small  scale  structure  with  a  wavelength  near  2  y 

in  the  2-7  km  layer.   This  structure  is  much  less  evident  in 

the  R.H.=602»  case  (Fig.  3.8a).    This  small  scale  structure 

is   near   the   region   of  maximum   rising   motion   above   the 

surface  front  position.   The  stronger  vertical  velocity  for 

the   R.H.=90%   case   (maximum   rising   motion   >  6xl0-4  cb/s 

versus  4xl0~4  cb/s  for  the  R.U. =60%    case)  is  consistent  with 

the  stronger  cumulus   convection.    This  convection   arises 

because  of  the  existence  of  conditional  instability  which  is 

due  to  the  prescribed  90%  relative  humidity.    The  cumulus 

convection   activity   appears   to   cause   the   small-scale 

variations  in  the  mixing  ratio  field. 

In   addition   to   a   larger   maximum,   the 

vertical  velocity  of  the  R.H.=90%  (Fig.   3.10b)  case  also 

displays  a  more   complex  structure,   with  two   concentrated 

regions  of  rising  motion.    The  primary  region  near  42°-43° 

extends  throughout  most  of  the  depth  of  the  troposphere  from 

about  1  km  to  11  km.    The  secondary  region  occurs  in  the 

middle-upper   troposphere   (5-11  km)   just   poleward   of   the 

primary   region   at   about   44°.     In   addition,   a   third, 
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relatively  weak,  maximum  upward  motion  occurs  farther 
poleward  near  47°.  This  structure  is  referred  to  as  a 
multi-band,  slanted  rising  motion  in  later  discussions.  This 
multi-band  structure  in  the  vertical  velocity  may  be 
compared  to  some  of  the  observed  multi-band  structure  of 
cloud  and  precipitation  associated  with  cold  fronts  during 
the  period  of  CYCLES  (CYCLonic  Ext  rat ropical  Storms)  Project 
(Parsons  and  Hobbs,  1983). 

As  may  be  expected,  the  stronger  warm  core 
and  stronger  convection  of  the  R.H . =  90%  case  is  related  to  a 
larger  latent  heat  release.  In  Fig.  3.14,  the  heating 
distribution  generally,  resembles  the  rising  motion  field 
(Fig.  3.10).  This  is  consistent  with  a  near  balance  between 
adiabatic  cooling  and  diabatic  heating.  The  R.H. =90%  case 
(Fig.  3.14b)  has  a  maximum  heating  rate  of  >6  °K/day,  which 
is  much  higher  than  the  1  °K/day  rate  of  the  R.H. =60%  case 
(Fig.  3.14a).  This  six-fold  increase  in  heating  is 
consistent  with  the  increase  in  precipitation  rate  shown  in 
Fig.  3.13.  The  greater  production  of  kinetic  energy  implied 
by  the  stronger  rising  motion  of  warm  air  gives  rise  to  a 
more  intense  horizontal  circulation  (Figs.  3.11b  versus 
3.11a).  The  upper  tropospheric  westerly  jet  (Fig.  3.11b) 
reaches  a  maximum  of  >  50  ra/s  (versus  the  >  35  m/s  for  the 
R.H. =60%  case)  near  52°,  which  is  two  degrees  poleward  of 
the  jet  core  of  the  R.H. =60%  case.   A  small  region  of  weak 

easterlies  that  develops  equatorward  of  the  westerly  jet  is 
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not  observed  in  the  R.H.=60%  case.  It  is  to  be  recalled 
here  that  these  are  perturbation  wind  fields;  the  basic  flow 
is  westerly  throughout  the  upper  troposphere.  The 
distribution  of  the  perturbation  zonal  wind  is  consistent 
with  an  anticyclonic  circulation  above  the  warm  core  that  is 
due  to  heating.  The  existence  of  this  easterly  regime 
limits  the  middle-upper  tropospheric  westerlies  to  poleward 
of  the  surface  front  position,  and  makes  the  middle 
tropospheric  westerly  regime  appear  as  a  narrow  jet-like 
structure  that  penetrates  downward  to  below  z  =  3  km.  The 
resultant  low  level  westerly  regime  is  sharper  and  more 
extensive  than  that  of  the  R.H.=60%  case. 

On  the  other  hand,  the  strength  of  the 
surface  easterlies  of  the  M60  and  M90  cases  is  very  similar. 
This  is  probably  because  most  of  the  impact  of  a  higher 
relative  humidity  is  in  deep  cumulus  convection,  with  the 
resultant  increased  latent  heat  release  occurring  mainly  in 
the  upper  and  middle  troposphere,  as  shown  in  Fig.  3.14. 

The  maximum  cyclonic  vorticity  is 
concentrated  within  the  frontal  zone  in  both  the  upper  and 
the  lower  troposphere  for  both  cases  (Fig.  3.12a  versus  Fig. 
3.12b).  The  stronger  horizontal  circulation  in  the  upper 
troposphere  for  the  R.H.=90%  case  is  also  evident  in  the 
vorticity  field. 
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(2)  Surface  Flux  Cases  (M61  versus  M91) 

When  surface   fluxes   of  momentum,   moisture 

and  heat  are  present,  the  differences  between  the  more  moist 

(R.H.=90%)   and   the   drier   (R.H.=60*)   cases   are   basically 

similar  to  those  between  the  no-flux  cases.    The  R.H.=90% 

case  has  a  warmer  upper  tropospheric  core,  as  represented  by 

the  decrease  of  height  of  the   0  =350°K  line,  and  again  has 

the  small-scale   variation   in   mixing  ratio   in   the   rising 

motion  region  over  the  surface  front  location   (Fig.   3 .  8d 

versus   Fig.   3.8c).   The   meridional   potential   temperature 

gradient  is  stronger  for  case  M91  than  for  case  M61  (Fig. 

3.9d   versus   Fig.   3.9c).     The   maximum   upward   motion, 

precipitation  and  heating  rate  are  much  larger  for  M91  than 

for   M61   (Figs.   3.10,    3.13   and   3.14).   Stronger   wind 

circulations  (Fig.  3. lid  versus  Fig.  3.11c)  are  shown  in  the 

R.H.=90%  case.     The   circulation   includes   a  more   intense 

westerly  regime  of  >  50  m  s_1  in  the  upper  troposphere,  the 

development  of  a  narrow  easterly  regime  equatorward  of  it, 

and  a  deeper  downward  extension  of  the  upper  westerlies  into 

the   low-level   easterly   jet.     The   differences   in   the 

magnitudes  of  vorticity  (Fig.  3.12d  versus  Fig.  3.12c)  are 

similar   to   other   fields   with   much   stronger   cyclonic 

vorticity  within  frontal  zone  for  the  R.H.=90%  case. 

It   is   to   be   noted   that   the   small-scale 

variations  in  all  fields  for  the  case  with  high  mixing  ratio 

and  surface  fluxes  (Fig.   3.8d-3.12d)  exist  near  surface  to 
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the  north  of  front.  As  can  be  seen  in  Fig.  3.15,  this"  is 
the  cold  air  region  where  the  sensible  and  latent  heat 
fluxes  are  large  and  apparently  cause  a  shallow  layer  of 
conditional  instability. 

On  the  other  hand,  there  is  one  noticeable 
difference  between  the  comparisons  of  the  with-flux  cases 
and  those  of  the  no-flux  cases.  The  multi-band  vertical 
velocity  structure  in  the  no-flux  (R.H.=90%)  case  (Fig. 
3.10b)  is  not  as  conspicuous  in  the  with-flux  (R.H.=90%) 
case  (Fig.  3.10d).  However,  a  weak  secondary  maximum  region 
does  appear  in  the  latter  near  7  km  at  46° ,  which  indicates 
a  double-maximum  structure  and  a  poleward  tilt  in  the 
vertical  structure.  The  physical  implication  of  this 
slantwise  structure  in  the  upward  motion  field  will  be 
discussed  later  in  Chapter  4. 

b.   Effect  of  Surface  Fluxes 

(1)  R.H.=60%  Cases  (M60  vs.  M61) 

The  inclusion  of  surface  fluxes  (Figs.  3.8a 

and  3.8c)  does  not  lead  to  any  discernible  differences  in 

potential  temperature  and  mixing  ratio  distributions,  except 

the  mixing  ratio  isolines  at  the  surface  have  been  spread 

towards   cold   air   when   surface   fluxes   are   included. 

Comparison   between   Figs.   3.9a   and   3.9c   shows   that   the 

surface  meridional  potential  temperature  gradient  for  the 

case  with  fluxes  substantially  weakens,  although  the  general 

pattern   delineating   the   frontal   zone   remains   the   same. 
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These  reductions  in  the  potential  temperature  and  mixing 
ratio  gradients  at  the  surface  are  due  to  the  distribution 
of  the  surface  fluxes  given  by  the  bulk  formula,  which 
provides  larger  heat  and  moisture  fluxes  poleward  of  the 
surface  front  in  the  cooler  and  drier  region  (Fig.  3.15). 

The  inclusion  of  surface  fluxes  causes  a 
slight  increase  in  the  vertical  velocity  (Fig.  3.10a  and 
3.10c).  This  is  expected  since  more  latent  heat  release 
occurs  with  the  increased  supply  of  heat  and  moisture  from 
below,  as  shown  in  the  comparison  between  Fig.  3.14c  and 
Fig.  3.14a.  Since  the  resultant  upper  tropospheric 
temperature  is  not  higher  (Figs.  3.8a  and  3.8c),  the 
additional  latent  heat  release  in  the  area  of  rising  motion 
is  compensated  almost  entirely  by  the  increased  adiabatic 
cooling  associated  with  the  stronger  upward  motion. 

Increased  conversion  of  available  potential 
energy,  as  implied  by  the  in-phase  relationship  between  the 
vertical  velocity  and  potential  temperature  fields,  leads  to 
a  stronger  upper  tropospheric  jet  (40  m/s  versus  35  m/s  in 
the  no-flux  case)  (Figs.  3.11a  and  3.11c).  On  the  other 
hand,  a  decrease  of  the  surface  easterlies  from  40  m/s  to 
30  m/s  occurs  due  primarily  to  the  dissipation  by  the 
momentum  flux  transfer.  The  maximum  easterly  jet  core  is 
raised  to  a  level  slightly  above  the  surface,  and  the 
strength   and   extent   of   the   narrow   lower   tropospheric 
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westerly  region  appear  to  be   slightly  reduced  for  no-flux 
case. 

Near  the  surface,  the  cyclonic  vorticity 
decreases  in  the  with-flux  cases  (Figs.  3.12a  and  3.12c). 
This  indicates  that  the  weakening  of  the  easterly  maximum 
and  thermal  gradient  of  the  front  at  the  surface  does  lead 
to  a  corresponding  reduction  in  the  horizontal  wind  shear  in 
the  lower  troposphere. 

(2)  R.H.=90%  Cases  (M90  vs.  M91) 

Similar  results  due  to  inclusion  of  the 
surface  fluxes  are  obtained  in  the  comparison  between  the 
with-flux  and  no-flux  cases  for  R .K. =90%    as  in  the  R.H.=60% 


cases 


1)  The  difference  in  the  distribution  of  potential 
temperature  and  mixing  ratio  is  not  discernible, 
except  near  the  surface  to  the  north  of  front  (Figs. 
3.8b  and  3. 8d) ; 

2)  At  the  surface,  the  inclusion  of  moisture  and  heat 
fluxes  weakens  the  meridional  gradients  of  mixing 
ratio  (Figs.  3.8b  and  3.8d)  and  potential 
temperature  in  the  frontal  zone  (Figs.  3.9b  and 
3.9d)  because  of  the  stronger  fluxes  in  the  cooler 
and  drier  regions; 

3)  The  surface  momentum  transfer  leads  to  a  weakened 
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surface  easterly  maximum  (Figs.  3.11b  and  3. lid), 
with  the  easterly  jet  core  occurring  slightly  above 
the  surface; 

4)  The  inclusion  of  surface  fluxes  causes  a  stronger 
and  more  concentrated  vertical  velocity  (Fig.  3.10) 
and  heating  rate  (Fig. 3. 14)  in  the  frontal  zone,  so 
that  the  upper  tropospheric  westerly  wind,  the 
westerly  regime  in  the  middle  and  lower  troposphere 
above  the  surface  front  are  stronger,  and  the 
meridional  wind  shear  in  the  lower  troposphere  is 
weaker  (Figs.  3.11b  and  3. lid  and  3.12b  and  3.12d); 
and 

5)  For  R.H.=90%,  the  multi-band  structure  in  vertical 
velocity  (Figs.  3.10b  and  3.10d),  for  the  with-flux 
case  is  less  conspicuous  than  that  for  the  no-flux 
case,  as  was  noted  earlier. 

C.   SUBTROPICAL  SIMULATIONS 

1.   Typical  Subtropical  Cases  (S90) 

To  examine  the  time  evolution  of  subtropical 
experiments,  the  case  of  S90  (initial  R.H.=90%  and  no 
surface  fluxes)  is  chosen  as  the  "typical  subtropical  case". 
The  initial  potential  temperature  and  mixing  ratio  are 
depicted  in  Fig.  3.16a  and  the  zonal  velocity  is  shown  in 
Fig.  3.16b.  The  choice  of  the  typical  subtropical  case  is 
motivated  by  an  intention  to  represent  the  major  differences 
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between  the  mean  midlatitude  atmosphere  and  the  mean 
subtropical  atmosphere  during  the  Mei-Yu  season.  The  mean 
subtropical  Mei-Yu  environment  is  conditionally  unstable 
with  a  R.H.  around  85%.  Furthermore,  during  intense 
convection  periods  when  a  LLJ  is  usually  observed,  the  R.H. 
is  around  90%  or  higher  (Chen  and  Tsay,  1978).  While  these 
intense  episodes  are  local  and  transient,  and  therefore 
cannot  be  properly  studied  by  the  present  two-dimensional 
model,  they  are  the  only  periods  when  the  LLJ  is  clearly 
defined.  Therefore,  the  somewhat  higher  value  of  R.H. =90% 
is  chosen  so  that  its  results  may  also  shed  light  to  the 
situations  when  an  intense  LLJ  is  developed.  In  this  view 
the  "typical  subtropical  case"  may  also  be  considered  to 
represent  the  "typical  intense  LLJ  case." 

As  in  the  typical  midlatitude  case,  the  large-scale 
deformation  forcing  causes  an  accumulation  of  the  potential 
temperature  and  moisture  isolines  to  result  in  a  frontal 
structure.  However,  the  approach  towards  a  quasi-steady 
state  is  slower.  Furthermore,  all  fields  develop 
significant  small-scale  variations  on  the  equator  side  of 
the  front.  (For  example,  see  the  potential  temperature  and 
mixing  ratio  fields  shown  in  Figs.  3.17a,b).  These 
variations,  which  are  similar  to  those  in  the  strong 
convection  cases  of  midlatitude  experiments,  are  apparently 
a  consequence  of  intense  cumulus  convection  in  the  model 
resulting  from  the  very  moist  and  unstable  atmosphere. 
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It  should  be  noted  that  once  the  initial  vertical 
0e  profile  is  set  from  the  temperature  sounding  and 
relative  humidity,  it  will  remain  fixed  in  time  at  lateral 
boundaries.  Since  a  conditionally  unstable  profile  occurs 
in  the  lower  latitude  regions  (because  of  the  initial 
potential  temperature  and  mixing  ratio  distribution),  the 
instability  will  persist  on  the  equatorward  boundary 
regardless  of  convective  adjustment  in  the  interior.  This 
slower  development  into  a  quasi-steady  state  is  probably  due 
to  two  factors:  the  weaker  large-scale  deformation  forcing 
(D=0 . 5xl0-5 s~ x )  prescribed  for  the  subtropical  cases  and  the 
presence  of  strong  conditional  instability  which  implies 
more  prominent  cumulus  convective  activity. 

The  small-scale  features  also  dominate  the  vertical 
velocity  field  such  that  the  frontal  structure  is  obscured 
by  high  frequency  oscillations.  To  remove  these  small-scale 
oscillations,  all  the  subtropical  experiments  fields  are 
subjected  to  a  daily  average  over  32  outputs  at  45  minute 
intervals.  Two  other  alternatives  were  considered:  to 
increase  the  control  of  large-scale  forcing,  or  to  decrease 
the  initial  relative  humidity.  However",  both  approaches 
will  reduce  the  contrasts  between  the  midlatitude  and  the 
subtropical  cases  especially  those  cases  with  clearly 
defined  LLJ,  which  is  one  of  the  main  objectives  of  this 
study.  Therefore,  it  was  decided  that  the  differences 
between  the  midlatitude  and  subtropical  conditions  should  be 
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kept  as  large  as  possible  to  identify  the  underlying 
dynamics  of  the  development  of  the  LLJ.  Intermediate  cases 
with  R.H.=70%  and  80%  will  also  be  run  to  ascertain  the 
physical  significance  of  the  intensively  convective  R.H.=90% 
cases . 

The  daily  averages  of  the  potential  temperature  and 
mixing  ratio,  and  the  meridional  potential  temperature 
gradient  for  days  2-5  are  shown  in  Figs.  3.18  and  3.19, 
respectively.  Here  each  "day"  represents  the  average  of  the 
previous  24  hours  (i.e.,  day  5  is  the  average  between  96- 
120  h)  .  In  these  diagrams,  the  small-scale  features  have 
been  filtered  out,  and  the  development  of  the  frontal 
structure  is  clearly  evident.  The  surface  mixing  ratio 
gradient  and  range  agree  reasonably  well  with  Chen's  (1978) 
observations  (Fig.  1.4b). 

The  meridional  potential  temperature  gradient  also 
evolves  into  a  frontal  structure  near  the  surface  around  16° 
after  4  days  (Fig.  3.19).  The  tilt  of  the  frontal  zone  is 
similar  to  the  midlatitude  case  (Fig.  3.4),  but  the  well- 
defined  frontal  zone  is  confined  to  the  lower  troposphere. 
Chen  and  Chang's  (1980)  comparison  between  different 
segments  of  an  elongated  Mei-Yu  front  that  extends  from  the 
subtropical  southern  China  to  the  midlatitude  central  Japan. 
They  observed  that  the  subtropical  segment  of  the  front  has 
less  vertical  tilt  and  is  shallower  than  the  midlatitude 
segment.   These  differences  were  attributed  to  the  stronger 
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convection  in  the  subtropics.  It  will  be  shown  later  that 
in  the  subtropical  with-flux  case  (S91)  the  tilt  is  indeed 
less  while  the  vertical  extent  remains  shallow. 

By  day  2,  the  concentrated  rising  motion  (Fig.  3.20) 
develops  in  the  equatorial  latitudes  and  then  moves  poleward 
towards  the  area  of  maximum  deformation  forcing.  At  day  3, 
the  maximum  rising  motion  is  at  16°  with  a  magnitude 
~17xl0~4  cb/sec.  After  day  4,  it  shifts  to  18°  and  the 
maximum  decreases  to  4-5  x  10~4  cb/sec.  This  is  the  process 
for  the  model  to  develop  the  Mei-Yu  front  from  an  initial 
zero-perturbation  state.  In  the  real  atmosphere,  the 
process  is  much  slower,  where  the  monsoon  trough,  which  is 
developed  over  a  seasonal  time  scale,  is  already  in  place 
when  the  front  moves  into  the  Mei-Yu  position  from  the 
north.  The  development  of  a  slanted  structure  in  the  lower 
(z=0-3  km)  and  upper  parts  (5-10  km)  of  the  concentrated 
rising  region  is  evident.  The  slanting  is  northward  with 
height,  in  the  same  sense  as  the  slope  of  the  frontal 
surface.  The  multi-band  structure  is  not  clearly  seen  in 
the  24-h  averaged  vertical  velocity  field  but  is  more 
clearly  visible  when  instantaneous  (not  daily-averaged) 
vertical  velocity  fields  are  examined.  An  example  is  given 
in  Fig.  3.21,  which  is  the  instantaneous  vertical  velocity 
at  the  time-step  of  day  4.5. 

The  evolution  of  the  zonal  wind  is  illustrated  in 
Fig.   3.22.    The  gross  structure  of  an  upper  tropospheric 
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westerly  jet,  poleward  of  a  lower  tropospheric  easterly  jet 
with  a  maximum  at  the  surface,  roughly  resembles  that  of  the 
typical  midlatitude  case  (M60,  Fig.  3.6).  However,  there 
are  significant  differences.  The  subtropical  westerly  jet 
core  is  lower  (at  10  km),  shallower  and  wider  in  the 
meridional  direction.  These  features  are  different  from 
what  would  be  expected  for  the  mean  polar  front  jet  and 
subtropical  jet.  Compared  to  Chen's  (1978)  observation  of  a 
subtropical  Mei-Yu  case  (Fig.  1.4a),  the  jet  core  is  lower. 
The  reason  for  the  lower  jet  core  (at  10  km  versus  12  km  in 
Fig.  1.4a)  is  apparently  the  fixed  tropopause  which  is 
imposed  by  the  mean  sounding  at  10  km.  In  an  intense  Mei-Yu 
period,  the  deep  convection  and  the  tropopause  are  a  couple 
of  kilometers  higher,  where  the  jet  maximum  is  also  located. 
On  the  other  hand,  the  meridional  extent  of  the  simulated 
westerly  regime  and  that  observed  are  comparable.  A 
significant  easterly  regime  occupies  the  entire  tropical 
upper  troposphere  and  caps  a  westerly  regime  in  the  tropical 
lower  troposphere.  In  Fig.  1.4a,  the  observed  zonal  wind  is 
also  easterly  in  the  upper  level  tropics,  although  the  area 
is  somewhat  higher  and  more  equatorward  than  the  simulation. 
The  imposed  tropopause  in  the  model  is  again  the  reason  for 
this  difference. 

Within  the  tropical  low-level  westerly  regime  is  a 
remarkably  strong  and  narrow  jet  structure  centered  at  z=3 

km  at  the  latitude  of  maximum  rising  motion.   This  low  level 
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jet  first  appears  at  day  2  with  a  magnitude  ~5  m/s.  It 
develops  to  >  10  m/s  at  day  5.  In  terms  of  the  relative 
position  with  respect  to  the  front,  this  low  level  jet  is 
similar  to  the  local  lower-tropospher ic  westerly  regime 
simulated  in  the  midlatitude  cases,  but  the  intensity  is 
much  stronger.  This  intensity  agrees  reasonably  well  with 
the  observed  LLJ  strength  in  Fig.  1.4b,  especially  when  the 
addition  of  a  mean  westerly  is  considered. 

The  strong  horizontal  wind  shear  that  develops 
around  the  frontal  zone  is  shown  in  Fig.  3.23.  Compared  to 
Fig.  3.7,  this  subtropical  case  has  a  deep  lower 
tropospheric  shear  zone  that  extends  into  the  middle 
troposphere,  and  its  magnitude  poleward  of  the  low  level  jet 
significantly  exceeds  that  of  the  upper  tropospheric 
maximum.  The  slope  of  the  frontal  shear  zone  is  again 
larger  (less  vertical  tilt)  in  the  subtropical  case  than  in 
the  midlatitude  case  (Fig.  3.7),  at  least  in  the  lower 
troposphere.  In  Fig.  3.23d,  a  region  of  strong  negative 
relative  vorticity  exists  above  the  surface  front  near 
10  km,  which  indicates  the  possibility  of  inertial 
instability.  This  is  the  consequence  of  the  intense 
convection  and  may  not  be  physically  realistic.  It  is  not 
known  whether  this  instability  contributes  to  the  small- 
scale  variations. 

To  demonstrate  that  the  five-day  integration  gives  a 
reasonable  quasi-steady  state,   the  averages  of  each  field 
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over  the  first  five  days  of  integration  was  used  as  the 
initial  conditions  for  another  integration  of  120  h.  The 
results  of  potential  temperature  and  mixing  ratio, 
meridional  potential  temperature  gradient,  vertical  velocity 
and  zonal  wind  for  the  (daily  averaged)  new  day  5  are  shown 
in  Fig.  3.24a-d.  These  fields  are  similar  to  those  of  the 
original  daily  averaged  at  day  5  (Figs.  3.18,  3.19,  3.20  and 
3.22).  This  result  indicates  that  the  daily  averaged  fields 
are  stable  in  time,  in  spite  of  the  presence  of  small-scale 
variations  and  possible  inertial  instability. 

The  subtropical  frontal  positions  as  indicated  by 
the  various  fields  are  farther  equatorward  from  the  central 
axis  of  deformation  forcing  than  that  of  the  midlatitude 
cases.  In  the  midlatitude  case,  the  central  axis  of 
deformation  is  at  46°  and  the  frontal  structure  near  the 
surface  develops  at  about  42°.  In  the  subtropical  case,  the 
axis  of  deformation  is  at  22.5°,  while  the  surface  frontal 
zone  is  near  16°.  This  increased  distance  apparently 
reflects  the  effect  of  the  smaller  Coriolis  parameter  which 
leads  to  a  relatively  larger  ageostrophic  secondary 
circulation  around  the  frontal  zone.  The  advection  by  the 
lower  branch  of  this  enhanced  secondary  circulation  tends  to 
move  the  surface  front  farther  equatorward.  In  addition, 
the  effect  of  a  possible  secondary  circulation  south  of  the 
front,  that  is  induced  by  the  strong  convection,  may  also 
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contribute  to  the  more  equatorward  position  of  the  surface 

front . 

2 .   Other  Subtropical  Cases 

As  in  the  midlat itude  experiments,   the  effects  of 

varying  the  initial  relative  humidity  and  surface  fluxes  are 

studied  by  running  three  additional  subtropical  cases  (S60, 

S61  and  S91).    The  resultant  structure  at  day  5,   averaged 

over  the  preceding  24-h  period,   for  all   four  subtropical 

cases  are  shown  in  Figs.  3.25-3.29.  The  precipitation  and 

heating  rates  are  shown  in  Figs.  3.30-3.31.   In  these  cases, 

no  significant  small-scale  variations  are  present,  so  the 

averaging  is  unnecessary.    Nevertheless,  it  is  done  to  make 

the  analysis  procedure  of  all  subtropical  cases  consistent. 

a.   Effects  of  Moisture 

Panels  a  and  b  in  Figs.   3.25-3.29  are  for  S60 

(R.H.=60%,   no-flux)   and   S90   (R.H.=90*.   no-flux)   cases, 

respectively.    The  R.H.=60%  case  (Fig.   3.25a)   has  a  less 

intense  moisture   front   throughout   the   lower   troposphere. 

The  poleward   vertical   tilt   of  the   maximum  mixing   ratio 

gradient  is  about  the  same  as  in  the  R.H.=90%  case  (Fig. 

3.25b).      However,   the   curvature   near   the   axis   of   the 

maximum  mixing  ratio  does  not  appear  as  sharp  as   in  the 

R.H.=90%  case  due  to  the  drier  air  and  weaker  gradient.   The 

short-  wave   feature  in   the  middle   troposphere  over   the 

surface  front  position  is  not  as  intense  as  in  the  R.H.=90% 

case,   but   it  is   still  very  notable,   which  indicates   the 
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presence  of  significant  cumulus  convection.  As  may  be 
expected,  the  R.H.=60%  case  has  a  lower  temperature  in  the 
upper  troposphere  (Figs.  3.25a  and  3.25b).  The  small-scale 
variations  near  the  surface  in  the  cold  air  (Fig.  3.25d- 
3.27d),  as  seen  in  midlatitude  cases,  evidence  the  shallow 
layer  of  conditional  instability  caused  by  surface  heat 
fluxes.  The  conditional  instability  can  be  examined  further 
by  calculating  the  9e  field  which  will  be  shown  later  in 
Fig.  4.6  (midlatitude  cases)  and  4.7  (subtropical  cases). 
The  weaker  meridional  potential  temperature  gradient  in  the 
frontal  zone  can  be  seen  by  comparing  Figs.  3.26a  and  3.26b, 
both  of  which  also  show  less  small-scale  features  elsewhere 
in  the  domain.  The  vertical  tilt  of  the  frontal  zone  is 
somewhat  larger  and  the  frontal  zone  extends  to  a  slightly 
lower  height  at  the  poleward  boundary. 

The  vertical  velocity  of  the  R.H.=60%  case  (Fig. 
3.27a)  has  a  concentrated  zone  of  rising  motion  in  the  lower 
and  middle  troposphere,  but  the  magnitude  is  less  than  in 
the  R.H.=90%  case. 

The  zonal  wind  of  the  R.H.=60%  case  is  much 
weaker  than  in  the  R.H.=90%  case  (Fig.  3.28a  versus  Fig. 
3.28b).  The  upper  tropospheric  westerly  jet  maximum  of 
<  25  m/s  is  only  about  70%  of  the  >  35  m/s  magnitude  of  the 
R.H.=90%  case.  This  is  consistent  with  the  expected 
difference  in  the  energy  conversion  between  the  two  cases, 
as  the  heating  rate  and  the  upward  motion  (which  correlates 
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almost  completely  with  the  heatin)g  of  the  warm  air  is 
stronger  for  the  R.H.=90%  case.  The  former  is  also  situated 
more  equatorward.  The  low-level  westerly  jet  structure 
still  persists,  but  the  magnitude  is  reduced  by  half.  The 
upper  tropospheric  easterly  regime  almost  disappears,  and 
the  surface  easterly  maximum  is  only  reduced  slightly  (Fig. 
3.28a  and  Fig.  3.28b).  These  features  also  appear  in  the 
comparison  between  Figs.  3.29a  and  b,  which  primarily 
indicates  a  significant  reduction  of  upper  tropospheric 
horizontal  wind  shear. 

When  fluxes  are  included,  the  comparison  between 
the  R.H.=60%  (S61)  and  R.H.=90%  cases  (S91)  [panels  c  and  d 
in  Figs.  3.25-3.31,  respectively]  yields  similar  results. 
In  general,  the  meridional  temperature  and  mixing  ratio 
gradients  and  the  zonal  wind  are  stronger  in  the  R.H.=90% 
"case. 

The  subtropical  R.H.=90%  cases  are  included  to 
study  the  intensely  convective  Mei-Yu  situations  when  the 
LLJ  is  fully  developed.  However,  the  high  initial  moisture 
and  the  subtropical  sounding  make  the  atmosphere  highly 
conditionally  unstable.  The  resultant  high  frequency,  small 
scale  variations  are  removed  by  the  24  h  averaging.  Even 
though  the  averages  remain  stable  for  extended  integrations, 
a  question  may  still  be  asked  as  to  whether  the  results  are 
not   distorted   seriously   by   the   small-scale   variations. 

Therefore,  two  intermediate  with-flux  cases,  R.H.=7  0%  (S71) 
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and  R.H. =80%  (S81),  are  run  to  assess  the  physical 
significance  of  the  S91  case.  Fig.  3.33  compares  the  24-h 
averaged  zonal  wind  of  all  four  subtropical,  with-flux 
cases.  Here  it  can  be  seen  that  the  pattern  of  all  four 
cases  are  similar.  The  midlatitude  upper  level  westerlies, 
lower  level  easterlies,  tropical  upper  easterlies  are 
present  in  all  cases,  and  the  intensity  of  these  features 
increases  systematically  as  the  initial  R.H.  increases  from 
60%  to  90%.  The  LLJ  is  indicated  by  a  very  small  area  of 
5  m/s  in  the  S61  case  (Fig.  3.33a),  this  area  expands 
considerably  with  increasing  humidity.  When  the  R.H. 
reaches  90%  (Fig.  3.33d),  the  jet  maximum  becomes  >  10  m/s. 
The  similarity  in  the  structure,  and  the  consistent 
intensification  of  the  circulation  with  increasing  moisture, 
indicate  that  the  subtropical  R.H. =90%  results  can  indeed  be 
viewed  as  physically  relevant  to  the  strong  LLJ  situations, 
b.   Effects  of  Surface  Fluxes 

When  surface  momentum,  heat  and  moisture  fluxes 
are  included  (Fig.  3.32),  these  important  contrasts  continue 
to  exist,  as  seen  by  comparing  panels  b  and  d  in  Figs.  3.25- 
3.31.  Similar  to  the  midlatitude  experiments,  the  addition 
of  surface  heat,  moisture  and  momentum  fluxes  leads  to  a 
higher  upper  tropospheric  temperature.  This  is  apparently 
due  to  the  increased  precipitation  and  the  associated  latent 
heat  release  (Figs.  3.30  and  3.31),  which  are  to  be  expected 
from  the  increased  sensible  and  latent  heat  supply   (Fig. 
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3.32).  A  weaker  potential  temperature  and  mixing  ratio 
gradient  at  the  surface  (Figs.  3.25  and  3.26)  is  due  to  the 
stronger  surface  fluxes  of  heat  and  moisture  in  the  cooler 
and  drier  region  (Figs.  3.32). 

The  vertical  velocity  of  the  with-flux  case 
displays  the  stacked  multi-band  structure  slanted  upward  and 
poleward  in  the  surface  frontal  zone  (Figs.  3.27b  and  d). 
Also,  the  magnitude  of  updraft  is  greater  than  in  the  no- 
flux  case.  This  is  consistent  with  the  precipitation  and 
heating  rates  (Figs.  3.30  and  3.31)  and  is  due  to  the 
additional  moisture  supply  from  the  lower  boundary  (Fig. 
3.32) . 

Comparisons  between  Figs.  3.28a  and  3.28c,  and 
between  Figs.  3.28b  and  3.28d  show  the  effects  of  surface 
fluxes  on  the  horizontal  circulation,  which  are  basically 
similar  to  the  results  of  the  raidlatitude  experiments.  When 
fluxes  are  included,  the  upper  tropospheric  zonal  winds, 
including  both  the  westerly  jet  and  the  less  intense 
easterly  jet,  are  enhanced.  The  lower  tropospheric  westerly 
jet  near  3-4  km  for  the  R.H.=90%  case  (Fig.  3.28d)  also 
intensifies.  Due  to  the  momentum  dissipation  at  the  surface 
(Fig.  3.32),  the  surface  easterlies  weaken,  and  the  jet  core 
is  shifted  slightly  upward  (Figs.  3.29b  and  d). 

The  comparisons  of  the  surface  fluxes  and  no 
fluxes  for  the  R.H.=60%  cases  (panels  a  and  c  of  Figs.  3.25- 
3.29)  basically  illustrate  similar  results. 
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D.   SUMMARY  OF  THE  RESULTS 

The  results  of  the  foregoing  numerical  experiments  may- 
be summarized  as  follows: 

1)  In  the  midlatitude  case,  the  zonal  velocity  field  is 
dominated  by  a  westerly  regime  occupying  almost  the 
entire  upper  troposphere  and  an  easterly  regime 
occupying  the  lower  troposphere.  On  the  other  hand,  in 
the  subtropical  cases,  the  zonal  flow  equatorward  of  the 
front  tends  to  reverse,  with  expansive  easterlies  in  the 
upper  troposphere  and  weak  westerlies  in  the  lower 
troposphere.  A  small  regime  of  easterlies  develops  in 
the  midlatitude  upper  troposphere  only  in  the  high- 
humidity  cases. 

2)  At  the  latitude  of  the  maximum  rising  motion,  a  narrow 
westerly  wind  regime  occurs  in  the  lower  troposphere. 
In  midlatitude  cases,  it  appears  as  a  small  downward 
(and  equatorward)  penetration  of  the  upper  tropospheric 
westerlies.  In  the  subtropical  cases,  it  exists 
distinctly  as  a  well-defined,  westerly  maximum  similar 
to  the  low-level  jet. 

3)  The  poleward  tilts  of  all  fields  except  vertical 
velocities  are  more  prominent  in  the  midlatitude  cases 
are   either   comparable    or   stronger    than   in    the 
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subtropical  cases.  The  frontal  zone  also  extends  higher 
in  the  midlat i tude  cases.  These  are  in  agreement  with 
observations  (e.g.,  Chen  and  Chang,  1980). 

4)  The  inclusion  of  higher  initial  humidity  values  or  of 
surface  fluxes  tends  to  give  stronger  vertical  velocity 
and  upper  tropospheric  zonal  winds.  This  is  a  result  of 
stronger  cumulus  convection  and  precipitation  whose 
latent  heat  release  provides  more  energy  input  to  the 
upper  troposphere.  At  lower  levels,  the  easteries  weaken 
as  a  consequence  of  the  momentum  dissipation  at  the 
surface . 

5)  When  the  vertical  velocity  is  intense,  a  multi-band 
structure  often  appears.  It  is  characterized  by  more 
than  one  concentrated  updraft  region,  stacked  closely  to 
each  other  along  the  poleward  slanted  frontal  surface. 
The  stronger  updraft  is  situated  equatorward  and  has  a 
larger  vertical  extent. 

6)  The  frontal  zone  develops  equatorward  of  the  dilation 

axis  of  the  large-scale  deformation  forcing.    In  the 

midlatitude  cases,  the  frontal  zone  is  closer  to  this 

axis   (about   4°   latitude   separation   at   the   surface), 

while  in  the  subtropical  cases  the  separation  is  greater 

(about  5-6°)  at  the  surface.  This  is  apparently  a  result 

of  the  smaller  f  and,   therefore,   a  relatively  larger 

ageostrophic   wind   which   advects   the   surface   front 

farther  equatorward.    In  addition,  the  effects  of  the 
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convection-induced   secondary    circulation   may    also 

contribute   to   the   more   equatorward   position   of   the 

surface  front. 

In  the  ensuing  section  these  numerical  results  will  be 
discussed  in  relationship  to  the  development  of  the  low- 
level  westerly  jet. 
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Fig.  3.1  Initial    potential    temperature     (solid),    0  (°K),    and    equivalent 

potential  temperature  (dash),  6e(°K),  soundings  at  the  center 
of  the  model  domain  for  midlatidude  cases  -  M60  and  M90 
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Fig.   3.15         Surface  fluxes  of  momentum   (dash;  -10-4    m  cb  s~2), 
sensible  heat   (solid;    -10~4    °K    cb   s"1)     and  latent 
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Fig.   3.16         As  in   Fig.   3.2,   except,  for   typical  subtropical  case,   S90. 
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Fig.  3.21         Cross  section  of  vertical  pressure  velocity,- a)(10"4    cb  s-1), 
at  day  4.5  for  case  S90. 
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Fig.   3.32  Surface  fluxes  of  momentum   (dash;   -10"4    m  cb   s^2), 

sensible  heat   (solid;      -10-4    °K    cb  s"1)     and   latent 
heat   (dot;  -10~4    <>K    cb  s"1).     (a)  S61  and   (b)   S91. 


114 


[\N*)Z 


si 


0»t    O'll    O'Ol     0"  8      0"  • 

i\n-A)z 


Of     o   I  oo 


CO 
o 

•1-H 

o 

Sh 
CO 


C 
CD 

t 

co 


o 


CO 

D 

03 

<D 
ft 
CO 

T3 

g     . 

•JH    | | 

_  CO 

co  ^ 

N 
CO 

bo 

cO 
u 

CO 


co 


00 


Bl  0>l  O'Ci  O'OI   OB    00     Ot     0  10   0 
(WX)Z 


cd  ^ 

0 

U  ^-* 

o 

3 

o 

o  _; 

111 

^ 

n 

CN  ^-. 

o 

3 

«m£ 

w 

h- 

O 

M 

1- 
< 

-1 

c  --1 

O  <£> 
23  W 

o 

o  _ 

oss-se 
ses   (a 

o 

U     CO 

„ 

o   0 

r» 

O 

CO 

o 

co 

p> 

CO 

be 

•i-i 

115 


IV.   MECHANISM  FOR  THE  LOW  LEVEL  JET  (LLJ) 

The  previous  section  shows  that  a  low-level  westerly 
regime  tends  to  develop  in  the  frontal  zone  due  to  the 
intensifying  meridional  temperature  gradient  in  the 
lowest  few  kilometers.  In  many  of  the  subtropical  cases 
this  low-level  westerly  regime  actually  develops  into  a 
low-level  jet  whose  height  of  2-4  km  is  in  good  agreement 
with  the  levels  (800-600  mb )  of  the  observed  low-level 
jet  (Chen,  1977;  1978).  The  jet  maximum  of  10-15  m/s  in 
these  subtropical  cases  is  also  in  reasonable  agreement 
with  observations. 

In  the  Introduction,  it  was  pointed  out  that  there  are 
several  theories  for  the  development  of  the  low-level  jet 
(LLJ).  Since  effects  of  terrain  (Holton,  1967;  Wexler, 
1961)  and  diurnal  cycles  (Blackadar,  1957)  are  not 
included  in  the  model,  the  related  boundary-layer 
theories  cannot  be  treated  in  this  study.  Neither  can 
the  theory  of  the  downward  transport  of  momentum  by 
cumulus  convection  (Masomoto  and  Ninomiya,  1969; 
Ninomiya,  1971),  because  the  convective  adjustment  used 
in  the  present  model  does  not  include  the 
parameterization  of  cumulus  momentum  transport. 

The  theory  proposed  by  Uccellini  and  Johnson  (1979)  in 
the  study  of  midlatitude  severe  storms  over  the  U.  S.  is 
based  on  a  wind-mass  adjustment,   whereby  the  secondary 
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circulation  forced  by  an  upper  level  jet  streak  causes  a 
vertical  coupling,  which  leads  to  the  development  of  a 
LLJ.  In  this  mechanism,  the  LLJ  will  develop  beneath  the 
upper  level  jet  streak  such  that  the  two  jets  will  be 
aligned  vertically  and  the  axes  will  be  perpendicular  to 
each  other.  Although  the  model-produced  low-level 
westerly  maximum  often  appears  (especially  in  the 
midlatitude  cases)  as  a  downward  extension  of  the  upper 
tropospheric  westerly  jet,  the  upper  level  jet  is  always 
situated  much  too  poleward  of  the  LLJ  core  for  wind-mass 
adjustment  mechanism  to  be  effective. 

In  the  present  study,  the  LLJ  in  the  subtropical 
experiments  is  much  better  defined  than  in  the 
midlatitude  experiments.  Furthermore,  in  all  cases  when 
the  vertical  motion  is  enhanced  due  to  increased  relative 
humidity  (therefore  also  enhanced  conditional 
instability)  and/or  surface  fluxes,  the  low-level 
westerlies  become  stronger.  This  is  also  true  with  the 
development  of  the  upper  level  easterlies  equatorward  of 
the  front,  which  also  appears  with  small  magnitudes  in 
the  midlatitude  experiments  with  higher  relative 
humidity.  In  all  subtropical  experiments,  these  upper- 
level  easterlies  expand  and  dominate  most  of  the  upper 
troposphere  equatorward  of  the  front.  These  results 
suggest  that  when  the  vertical  motion  above  frontal  zone 
is  intense,  a  "reversed  Hadley  cell"  exists  equatorward 
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of  the  front,  with  rising  motion  in  the  frontal  region 
and  sinking  motion  further  equatorward.  Such  an 
overturning  would  provide  an  easterly  Coriolis  torque  at 
the  upper-level  equatorward  flow  and  a  westerly  Coriolis 
torque  at  the  lower-level  poleward  flow,  which  may 
explain  the  upper-level  easterlies  and  the  low-level 
westerly  jet.  Indirect  evidence  of  the  "reversed  Hadley 
cell"  to  the  south  of  LLJ  may  be  found  in  an 
observational  study  by  Ninomiya  and  Akiyama  (1974).  They 
found  that  moist  air  is  located  above  the  LLJ  in  the 
middle  and  upper  troposphere,  but  the  air  south  of  the 
LLJ  is  fairly  dry  throughout  the  whole  troposphere.  This 
is  indicative  of  a  large-scale  subsidence  south  of  the 
front . 

More  evidence  of  this  "reversed  Hadley  cell"  is  seen 
in  the  temperature  distribution  of  all  cases  in  which 
strong  cumulus  convection  is  indicated.  For  example,  the 
temperature  in  Fig.  3.18  (S90  -  the  typical  subtropical 
case)  increases  equatorward  throughout  the  troposphere 
for  the  first  day.  Afterward,  reversed  meridional 
temperature  gradients  begin  to  appear  in  the  tropical 
region,  such  that  by  day  5  the  maximum  temperature  at  any 
level  in  the  middle-upper  troposphere  is  near  the  frontal 
zone.  This  situation  is  also  apparent  in  Fig.  3.25  for 
all  subtropical  cases.  The  reversed  temperature  gradient 
to   the   south   of   the   maximum   upward   motion   has   been 
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observed  both  in  the  Mei-Yu  front  (Chen  and  Chang,  1980) 
and  in  the  Baiu  front  (Matsumoto  et  al.,  1970;  Ninomiya 
and  Akiyama,  1970).  Therefore,  the  equatorside  "reserved 
Hadley  cell"  is  actually  thermally  direct,  with  warm  air 
rising  in  the  frontal  zone  and  cool  air  sinking 
equatorward.  Hence,  this  circulation  can  be  maintained 
and  as  well  as  continue  to  provide  the  upper-level 
easterly  and  lower-level  westerly  accelerations. 

The  next  question  concerns  the  identification  of  the 
main  mechanism(s)  for  this  equatorward  overturning. 
There  appears  to  be  three  possible  mechanisms: 

1)  The  secondary  circulation  normally  associated  with 
the  large-scale  flow  of  a  front; 

2)  The  cumulus  convection  due  to  the  conditional 
instability,  which  can  increase  the  latent  heat 
release  and  force  an  enhanced  secondary 
circulat  ion ; 

3)  The  third  possibility  is  based  on  the  observation 
that  the  vertical  velocity  field  of  several  of  the 
more  intense  cases  exhibits  a  multi-band 
structure.  Each  of  the  rising  branches  is 
upright,  and  is  stacked  in  a  poleward  slanting 
configuration.  In  a  way  this  structure  resembles 
the  "slantwise  convection"  associated  with 
symmetric  instability  in   a  baroclinic   atmosphere 
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(e.g.,  Emanuel,  1985).  This  slantwise  convection 
may  also  be  a  relevant  mechanism. 

To  differentiate  the  last  two  mechanisms,  the  term 
"upright  convection"  will  be  used  to  describe  the 
basically  single-band,  vertically-oriented  convection 
which  is  normally  associated  with  conditional 
instability,  and  "slantwise  convection"  will  be  used  to 
describe  the  slanted  and  stacked  multi-band  structure 
which  resembles  the  convection  predicted  by  the  symmetric 
instability  theory.  Fig.  4.1  illustrates  schematically 
the  differences  among  the  three  mechanisms. 

The  meridional-vertical  s treamf unct ion  is  shown  in 
Figs.  4.2  and  4.3,  for  the  midlatitude  and  subtropical 
cases,  respectively.  For  the  midlatitude  R.H.=60%  cases 
(M60  and  M61),  the  atmosphere  remains  statically  stable, 
and  Figs.  4.2a,c  indicate  that  the  large-scale  frontal 
secondary  circulation  consists  of  a  single  cell  in  the 
Hadley  sense.  When  R.H.=90%,  the  initial  sounding  of  the 
midlatitude  cases  (M90  and  M91)  is  conditionally  unstable 
in  the  lower  latitudes.  At  day  4,  both  cases  show  a 
small  and  very  weak  "reversed  Hadley  cell"  in  the  middle 
troposphere  equatorward  of  the  front.  Apparently  this 
cell  is  responsible  for  the  small  region  of  upper  level 
easterlies  equatorward  of  the  front.  However,  the  modest 
cumulus   convection   through   the   release   of   conditional 
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instability  is  not  intense  enough  to  force  a  strong" 
reversed  Hadley  ceil  that  can  develop  a  strong'  LLJ. 

All  the  subtropical  cases  are  initially  conditionally 
unstable.  A  larger  and  stronger  "reversed  Hadley  cell" 
is  present  in  all  cases  (Fig.  4.3).  This  cell  is 
situated  poleward  and  is  more  vertically  confined  for 
the  R.H.=60%  case  than  for  the  R.H.=90%  case  with  the 
middle-level  rising  motion  occurring  poleward  of  the 
surface  front.  This  is  the  situation  when  the  LLJ  is 
either  ill-defined  (S60  case,  Fig.  3.28a)  or  rather  weak 
(S61,  Fig.  28c).  When  R.H.=90*,  the  reversed  Hadley  cell 
is  located  more  equatorward,  is  more  intense  and  extends 
throughout  the  vertical  domain  (Figs.  4.3b,d).  The  LLJ 
is  well-defined  in  this  situation  (Figs.  3.28b,d),  and 
clearly  coincides  with  the  concentrated  poleward  flow  at 
the  lower  branch  of  the  reversed  Hadley  cell.  The 
cumulus  convection  is  the  strongest  in  these  cases  and  is 
apparently  responsible  for  generating  the  intense 
equatorside  meridional  circulation  which,  in  turn,  leads 
to  the  development  of  the  LLJ. 

It  is  clear  that  the  cross-frontal  secondary 
circulation  is  different  from  the  secondary  circulation 
equatorward  of  the  front,  and  therefore  it  can  not  induce 
the  LLJ.  However,  the  roles  of  upright  convection  and 
slantwise  convection  are  not  easy  to  distinguish.  The 
slantwise   convection    structure   is    present   in   both 


121 


midlat i tude  (when  the  LLJ  is  not  present)  and  subtropical 
(when  the  LLJ  is  developed)  cases.  The  slanting  is 
apparently  related  to  the  structure  of  the  frontal 
surface.  Although  the  frontal  forcing  cannot,  by  itself, 
directly  cause  the  equatorward  secondary  circulation,  it 
exerts  a  strong  influence  on  the  environmental  conditions 
that  affect  convection.  This  may  involve  two  mechanisms. 
The  first  is  the  large-scale  lifting  of  moist  air.  This 
leads  to  either  stable  precipitation  or  a  "CISK" 
(Conditional  Instability  of  the  Second  Kind)  type 
convection,  where  the  low-level  convergence  is  provided 
by  the  frontal  forcing.  The  second  involves  the  creation 
of  conditions  for  symmetric  instability.  These  conditions 
occur  when  the  vertical  gradient  of  equivalent  potential 
temperature  (  8e)  becomes  negative  along  the  surfaces  of 
pseudo-angular  momentum,  which  is  defined  by  the 
following  function  (Emanuel,  1983a, b): 

M  =  u  +  fy. 

In  a  symmetrically  unstable  baroclinic  atmosphere, 
slantwise  convection  occurs  when  there  is  a  slantwise 
displacement  of  an  air  parcel  along  a  constant  M  surface. 
Based  on  a  diagnostic  semi-geos trophic  model,  Emanuel 
(1985)  further  showed  that  even  in  the  presence  of  small 
symmetric  stability,  frontal  forcing  can  result  in  a 
strong,   concentrated,   sloping  updraft  which  is   located 
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ahead  of  the  region  of  maximum  geostrophic  fronto- 
genetical  forcing.  This  possibility  is  supported  by 
Sanders  and  Bosart  (1985)  who  reported  observations  of  a 
concentrated  sloping  updraft  that  resembles  slant-wise 
convection  in  a  mesoscale  precipitation  band.  The  sloping 
updraft  remains  even  after  the  environment  becomes 
symmetrically  stablized.  Apparently  the  f ronto-genet ical 
process  can  force  the  moist  updraft  along  the  sloped  M- 
surfaces  where  narrow  convective  bands  can  develop. 
Thus,  slantwise  convection  may  be  expected  in  areas  of 
neutral  or  small  symmetric  stability  when  f rontogenet ical 
forcing  is  present. 

The  evolution  of  6e  and  the  M-function  for  the 
typical  midlatitude  and  subtropical  cases  (M60  and  S90) 
is  shown  in  Figs.  4.4  and  4.5.  In  the  midlatitude  case, 
the  f rontogenet ical  area  is  statically  stable  initially 
and  the  M-surfaces  tilt  only  slightly  (Fig.  4.4a).  As 
the  f rontogenet ical  process  proceeds,  the  frontal 
structure  develops  and  leads  to  a  stronger  tilt  of  the  M- 
surfaces  in  the  highly-barocl inic  frontal  zone.  This 
results  in  an  area  of  near-neutral  symmetric  instability 
after  day  1,  although  the  area  is  limited  to  the  lowest 
troposphere.  In  the  subtropical  case,  conditional 
instability  is  present  throughout  the  domain  initially 
(Fig.  4.5a).  This  condition  continues  until  day  4  when 
the  frontal  structure   is  developed  (Fig.   3.19)   and  the 
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frontal  zone  becomes  stable  in  the  upright  convection 
sense.  As  a  result  of  f rontogenes is ,  the  M-surfaces  tilt 
more  towards  the  cold  air  in  the  frontal  zone.  This 
increased  tilt  leads  to  a  condition  of  symmetric 
instability  in  a  relatively  deep  layer  between  the 
surface  and  6  km  in  the  frontal  zone.  As  the  convective 
instability  is  eliminated,  a  slantwise  convection 
structure  develops  in  the  vertical  motion  field.  It 
appears  that  this  structure  is  produced  by  the 
front ogenesis  forcing  of  moist  updrafts  along  the  sloped 
M-surfaces,  similiar  to  Sanders  and  Bosart's  (1985) 
observat i  ons . 

Distribution  of  day-5  9e  and  the  M-function  are 
presented  in  Figs.  4.6  and  4.7  for  the  midlatitude  and 
subtropical  experiments,  respectively.  In  the  R.H.=60% 
midlatitude  cases  (M60  and  M61,  Figs.  4.6a,c),  static 
stability  remains  positive  and  the  area  favorable  for 
slantwise  convection  is  rather  limited  and  is  confined  to 
the  lowest  2  km  in  the  frontal  zone.  Thus,  the  secondary 
circulation  is  almost  entirely  a  result  of  the  frontal 
forcing. 

Distribution  of  0e  and  the  M-function  are  given  in 
Figs.  4.6b,d  for  the  R.H.=90%  midlatitude  cases  (M90  and 
M91).  In  these  cases,  the  equatorward  side  of  the  front 
is  conditionally  unstable  while  the  poleward  side  is 
stable  initially.    It  is  interesting  to  note  that  in  the 
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frontal  zone  the  convective  adjustment  apparently  has 
stablized  the  atmosphere  in  the  vertical  direction. 
However,  the  0  ?  and  M  lines  are  nearly  parallel  in  the 
frontal  zone,  an'd  the  area  favorable  for  sloping  updraft 
is  much  greater  than  that  for  the  R.H.=60%  cases.  In 
this  case,  a  low  level  maximum  wind  of  5  m/s  actually 
develops  near  z=3  km. 

Similar  to  the  R.H.=90%  midlatitude  cases,  all  the 
subtropical  without-flux  cases  (Fig.  4.7a,b)  have 
conditional  unstable  profiles  equatorward  of  the  front 
and  a  stable  9e  profile  poleward.  In  the  with-flux  cases 
(Fig.  4.7c,d),  the  conditional  instability  also  exists  in 
the  poleward  region,  a  result  of  the  supply  of  moisture 
and  heat  at  the  surface.  The  frontal  rising  motion  zones 
in  all  cases  again  have  stable  6  e  profiles,  but  shows 
signs  of  even  stronger  symmetric  instability  in  lower 
troposphere . 

In  the  S60  case  (Fig.  4.7a),  there  is  a  larger  area  of 
near  neutral  9e  profiles  along  the  M-surfaces  in  the 
vicinity  of  the  frontal  rising  motion.  In  the  S61  case 
(Fig.  4.7c),  the  9  e  profiles  in  the  lower  troposphere 
actually  become  unstable  along  the  M  surfaces  in  the 
rising  motion  zone.  The  same  is  true  for  the  S90  and  S91 
cases  (Figs.  4.7b,d),  except  the  area  of  symmetric 
instability  is  even  larger  than  in  the  R.H.=60%  cases. 
Thus,   for   all   cases   when   the   condition   of   symmetric 
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instability  is  suggested  by  the  9  e  profile  along  M- 
surfaces,  a  relatively  intense  LLJ  is  developed.  While 
in  the  present  model  framework,  it  is  difficult  to  assess 
directly  the  possible  roles  of  symmetric  instability,  the 
appearance  of  the  slantwise  convection  structure  seems  to 
be  an  important  factor  in  the  occurrence  of  the 
equatorward  secondary  circulation. 
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Fig.   4.1  Schematic       diagram      to       show       the       differences       between 

(a)    frontal     secondary     circulation,     (b)     upright     convection, 
and   (c)   slantwise  convection. 
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Fig.   4.4  Cross    section    of    M-f unction     (solid)    in    unit    of    m    s-1      and 

equivalent    potential    temperature    (dash),     6e(°K)      for    typical 
midlatitude  case,  M60,  at   (a)   day  0.0,   (b)   day   1.0. 
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Fig.  4.5  Cross    section    of    24-hour    averaged     M- function     (solid)     in 

units   of  m   s_1     and   equivalent  potential   temperature    (dash), 


Se(°K)      for     typical 
(b)    day   1.0. 


subtropical     case,    S90,     at     (a)     day     0., 
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V.   SUMMARY  AND  CONCLUSIONS 

In  this  work  a  two-dimensional  frontal  model  is  used  to 
study  the  structure  and  behavior  of  the  Mei-Yu  front,  with 
special  emphasis  on  the  dynamics  of  the  low-level  jet  (LLJ) 
that  is  frequently  observed  ahead  of  the  front.  The  Mei-Yu 
front  develops  in  the  area  of  the  east  Asian  monsoon  trough 
during  early  summer,  and  is  the  most  important  rain- 
producing  system  over  east  Asia  outside  of  the  typhoon 
season.  This  study  is  motivated  by  the  unique  mixed 
midlat i tude-t ropical  properties  of  the  Mei-Yu  front  (e.g., 
Chen  and  Chang,  1980),  where  heavy  convection  and  strong 
low-level  horizontal  shear  are  associated  with  the  nearly 
zonally-oriented  front.  The  occurrence  of  the  LLJ  is  highly 
correlated  with  heavy  convective  rainfall.  The  mechanism  for 
the  development  of  this  LLJ  is  the  main  focus  of  this  study. 

In  the  numerical  model,  the  time-varying  response  to  a 
constant  large-scale  forcing  was  treated  by  a  set  of 
perturbation  equations.  The  large-scale  forcing  was 
prescribed  by  a  stretching  deformation  field  which  is 
characteristic  of  the  f rontogenet ical  forcing  of  the  Mei-Yu 
front.  The  quasi-steady  state  responses  to  this  forcing 
were  obtained  by  integrating  the  perturbation  equations  from 
an  initial  state  which  is  given  by  the  seasonally-averaged 
zonal  flow.  Two  major  sets  of  experiments  were  conducted. 
The   first   is   based   on  midlatitude   conditions,   where   the 
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Coriolis  parameter  at  46°   latitude  and  an  initially  stable 

sounding  were  chosen.    The  second  is  based  on  subtropical 

conditions,  where  the  Coriolis  parameter  at  22.5°   latitude 

and  an  initial  sounding  that  is  conditionally  unstable  was 

chosen.    For  each  set  of  initial  conditions,   a  number  of 

supplemental  cases  were  integrated  to  study  the  effects  of 

varying   initial   humidity   values   and   surface   fluxes   of 

momentum,  heat  and  moisture.   In  the  midlatitude  cases,  the 

quasi-steady   state    is    reached   within    four    days    of 

integration.   In  the  subtropical  cases,  small-scale  features 

due   to   intense   convection   cause   the   motion   field   to 

fluctuate  throughout  the  integration.    A  24-hour  averaging 

was  applied  at  day  5  (96-120  h).    This  is  considered  the 

quasi-steady  state  as   further  integration   shows  that   the 

motion  remains  quasi-steady. 

For  both  midlatitude   and  subtropical  experiments,   the 

frontal  zone  develops  equatorward  of  the  dilation  axis.  The 

separation   between   the   front   and   the   dilation   axis   is 

greater  for  the  subtropical  cases  than  for  the  midlatitude 

cases.    This  larger   separation  is  apparently  due  to   the 

smaller  f  and,   therefore,  a  relatively  large  ageostrophic 

secondary  circulation  around  the  frontal  zone  whose  lower 

branch  advects  the  surface  front   farther  equatorward.   The 

midlatitude  front  extends  deeply  into  the  upper  troposphere 

with  a  strong  poleward  tilt,  whereas  the  subtropical  front 

is  confined  to  the  lower  troposphere  with  less  tilt,  which 
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is  in  good  agreement  with  observations  (e.g.,  Chen  and 
Chang,  1980).  Along  the  sloping  front,  slantwise  updrafts 
develop  with  a  multi-band  structure.  This  updraft  is  more 
evident  in  the  subtropical  cases  and  in  the  more  moist 
midlatitude  cases. 

A  westerly  jet  in  the  upper  troposphere  and  an  easterly 
jet  in  the  lower  troposphere  develop  in  both  midlatitude  and 
subtropical  cases,  with  smaller  intensities  in  the 
subtropical  cases.  The  inclusion  of  higher  humidity  or 
surface  fluxes  leads  to  stronger  upper  westerlies, 
apparently  as  a  result  of  stronger  cumulus  convection.  The 
lower  level  easterlies  weaken  both  in  the  midlatitudes  and 
the  subtropics  with  the  inclusion  of  surface  fluxes.  This 
is  apparently  due  to  the  effect  of  momentum  dissipation  at 
the  surface.  An  interesting  feature  in  the  zonal  wind  field 
in  the  subtropical  cases  is  the  concurrent  development  of 
upper-level  easterlies  and  low-level  westerlies  equatorward 
of  the  front.  The  low-level  westerly  maximum  at  z  =  3-4  km 
resembles  a  LLJ,  whose  intensity  increases  when  higher 
moisture  is  included.  In  the  midlatitude  cases,  only  a 
small  downward  and  equatorward  penetration  of  the  upper 
westerlies  appears  in  the  similar  area. 

While  several  theories  have  been  proposed  for  the 
development  of  low-level  jets  in  the  Mei-Yu  and  other 
circulation  systems,  the  choice  of  the  model  set  up  prevents 
consideration   of   some   of   these   theories   as   a   possible 
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explanation  for  the  simulated  LLJ.  The  absence  of  terrain 
and  diurnal  radiation  cycles  in  the  model  precludes  the 
boundary-layer  theories  (Wexler,  1961;  Holton,  1967; 
Blackadar,  1957),  and  the  absence  of  a  cumulus  momentum 
transport  parameterization  precludes  the  mechanism  due  to 
cumulus  momentum  transport  (Masomoto  and  Ninomiya,  1969; 
Ninomiya,  1971).  In  addition,  the  vertical  coupling 
mechanism  suggested  by  Uccellini  and  Johnson  (1979)  also 
does  not  apply  in  the  Mei-Yu  simulation,  because  the  upper 
tropospheric  westerly  jet  is  located  too  far  from  the 
latitude  of  the  LLJ  (Fig.  1.4). 

The  appearance  of  an  easterly  regime  in  the  upper 
troposphere  on  the  equator  side  of  the  front  above  the  LLJ 
suggests  that  the  LLJ  may  be  the  result  of  a  secondary 
circulation  that  resembles  a  "reversed  Hadley"  cell.  This 
possibility  is  supported  by  the  temperature  distribution, 
which  indicates  that  the  "reversed  Hadley"  cell  is  actually 
a  thermally  direct  circulation  that  is  capable  of 
maintaining  the  kinetic  energy.  The  Lower  branch  of  this 
circulation  is  directed  poleward  toward  the  front,  and  it 
exerts  a  westerly  Coriolis  torque  which  can  develop  and 
maintain  the  LLJ.  A  meridional-vertical  streamf unction 
reveals  the  structure  of  this  secondary  circulation,  with  a 
strong  lower  branch  return  flow  coinciding  with  the 
development  of  a  LLJ  in  the  more  moist,  subtropical  cases. 
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Several  possible  mechanisms  may  generate  the  reversed 
Hadley  cell.  The  possibility  that  the  reversed  Hadley  cell 
is  a  secondary  circulation  associated  with  the  basic  frontal 
structure  was  ruled  out  because  in  the  drier  midlatitude 
cases  the  frontal  secondary  circulation  is  a  single  cell 
around  the  front  with  upward  motion  on  the  equator  side  and 
downward  motion  on  the  pole  side.  No  LLJ  is  developed  in 
these  cases.  In  the  more  moist  (especially  the  subtropical) 
cases,  the  LLJ  developed  with  the  equator-side  secondary 
circulation.  This  suggests  the  importance  of  cumulus 
convection  in  developing  the  reversed  Hadley  cell  and  the 
LLJ. 

In   all   cases   where   a   LLJ   is   developed,   there   is   a 

sloping  updraft  with  multi-band  rising  motion  cells.    This 

suggests  that,   besides  regular  "upright"  convection  due  to 

conditional  instability,  slantwise  convection  arising  out  of 

symmetric  instability  may  also  be  relevant.   The  roles  of 

these  two  types  of  convection  were  studied  by  examining  the 

evolution    of    the    equivalent    potential    temperature 

distribution  in  the  vertical  and  along  the  pseudo-angular 

momentum  (M-function)   surfaces.    In  the  subtropical  cases 

where  the  atmosphere  is  conditionally  unstable  initially, 

upright   convection   develops   rapidly   and   the   conditional 

instability  is  quickly  released.   During  the  early  stages  of 

the  development,  the  M-function  surface  is  nearly  vertical, 

so  that  the  condition  for  symmetric  instability,  (99e/8z)M  < 

140 


0,  is  also  given  by  that  for  conditional  instability,  (99©/ 
3z)y  <  0.  As  the  front  ogenet  ical  process  continues,  the  M- 
surfaces  in  the  frontal  zone  begin  to  slope  and  the 
condition  for  symmetric  instability  become  different  from 
that  for  conditional  instability.  As  the  conditionally 
unstable  soundings  adjusts  to  neutral,  the  0e  sounding 
along  M  surfaces  is  either  neutral  or  slightly  unstable  in 
the  middle-lower  troposphere  below  6  km.  At  this  stage, 
slantwise  convection  is  observed  in  the  frontal  zone  and  is 
indicative  of  the  possible  importance  of  symmetric 
instability  processes. 

In  the  more  moist  midlatitude  cases,  slantwise 
convection  is  also  observed.  In  these  cases,  the  frontal 
zone  remains  statically  stable.  However,  the  f rontogenet ic 
forcing  produces  slanted  M  surfaces  there,  with  the  0e 
profiles  on  these  surfaces  reduced  to  only  slightly  stable 
based  on  the  symmetric  instability  criterion.  This  result  is 
in  agreement  with  the  recent  theory  of  Emanuel  (1985)  who 
showed  that,  under  f rontogenet ic  forcing,  sloping  updrafts 
can  occur  in  a  concentrated  region  on  the  warm  side  of  the 
frontal  zone,  and  result  in  an  enhanced  secondary 
circulation.  It  also  agrees  with  the  observation  by  Sanders 
and  Bosart  (1985)  of  a  mesoscale  precipitation  band  in  a 
winter-time  extratropical  cyclone  over  northeastern  United 
States,   in   which   a   sloping  updraft   was   observed   in   an 
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environment   in   which   the   symmetric   stability  was   weakly 
pos  i t  ive . 

The  main  findings  of  this  study  are  that  the  LLJ  in  a 
subtropical  Mei-Yu  front  can  develop  as  a  result  of  the 
Corioiis  torque  of  a  secondary  circulation  on  the  equator 
side  of  the  front.  This  secondary  circulation  is  thermally 
direct  and  is  in  the  reversed  Hadley  sense.  It  is  driven  by 
intense  condensation  heating  due  to  heavy  precipitation. 
The  upward  motion  takes  the  form  of  slantwise  convection  as 
a  result  of  the  f rontogenes is  forcing.  The  slantwise 
convection  begins  to  appear  after  the  atmosphere  is 
stablized  by  the  convective  adjustment.  It  is  conceivable 
that,  without  the  slantwise  convection,  which  is 
concentrated  along  the  frontal  surface,  the  latent  heat 
release  would  spread  out  over  a  broader  area.  This  would 
give  rise  to  a  weaker  secondary  circulation.  These 
conclusions  are  consistent  with  the  observations  that  low- 
level  jets  are  often  associated  with  intense  convection  and 
heavy  rainfall  in  the  Mei-Yu  front,  and  the  indirect 
evidence  of  a  sinking  region  south  of  the  Baiu  front  as 
revealed  by  Matsumoto's  (1972)  moisture  analysis.  Further 
study  is  needed  to  clarify  the  possible  roles  of  the 
symmetric  instability  mechanism.  However,  the  strong  in- 
phase  correlation  between  the  rising  velocity  and  heating 
indicates  that  the  main  conclusion  concerning  the  role  of 
the  secondary  circulation  should  not  depend  on  the  details 

142 


of  the  formulation  of  condensation  heating  in  the  model.  In 
fact,  the  frontal  forcing  is  the  dominating  factor  that 
determines  the  region  of  maximum  rising  motion.  Given  the 
moist  subtropical  environment  for  the  Mei-Yu  system,  if  a 
specified  (instead  of  predicted)  heating  function  is 
appropriately  prescribed  in  the  model,  a  strong  secondary 
circulation  equatorward  of  the  front  may  also  be  expected. 

Observational  studies  (e.g.,  Chen  and  Tsay,  1978) 
reported  that  the  Mei-Yu  convection  and  heavy  rainfall  are 
organized  in  mesoscale  convective  complexes  (MCC),  which  can 
not  be  treated  in  the  present  two-dimensional  model.  A 
three-dimensional  mesoscale  numerical  predication  model, 
where  the  MCC ' s  can  be  explicitly  simulated,  should  be  used 
to  reexamine  the  relationship  between  the  LLJ,  secondary 
circulation  and  convection.  A  three-dimensional  model  is 
also  necessary  to  study  the  entire  life  cycle  of  the  Mei-Yu 
front  and  the  LLJ.  The  development  and  decaying  stages  of 
the  LLJ  and  the  Mei-Yu  front  are  of  different  time  scales. 
While  the  Mei-Yu  front  usually  persists  for  several  days  to 
one  week  or  more,  the  appearance  of  the  LLJ  is  more 
transient  in  nature.  The  time  evolution  of  the  LLJ  within 
the  quasi-steady  Mei-Yu  front  is  therefore  an  interesting 
topic  for  further  studies. 

The   model   may   also   be   improved   to   study   additional 

effects.   The   inclusion   of   a   planetary   boundary   layer 

parameterization  that   is  more  sophisticated  than  the  bulk 
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aerodynamics  and  is  coupled  with  a  diurnal  radiation  cycle 
will  better  represent  the  effects  of  friction,  especially 
the  interaction  with  radiation  processes.  The  inclusion  of 
terrain  may  be  particularly  important,  since  the  Mei-Yu 
front  is  over  the  mountainous  region  of  southern  China, 
Taiwan  and  Japan.  The  topographical  effect  of  the  mountains 
and  the  land-sea  complex  in  this  region  undoubtedly  play  an 
important  role  in  both  the  large-scale  front  and  the 
embedded  mesoscale  convective  systems.  In  addition,  Keyser 
and  Anthes  (1982,  1986)  have  argued  that  a  stratification- 
dependent  boundary  layer  effect  is  important  for  the 
intensity  of  vertical  velocity  in  a  surface  front,  while  Mak 
and  Bannon  (1984)  argued  that  condensation  heating  alone 
plays  a  more  critical  role.  A  sophisticated  planetary 
boundary  layer  parameterization  may  help  to  clarify  this 
issue.  Finally,  further  observational  studies  are  also 
necessary  to  investigate  the  Mei-Yu  processes  both  in  the 
planetary  boundary  layer  as  well  as  in  the  entire 
troposphere . 
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Appendix 

A.   SPECIFICATION  OF  MODEL  GRIDS 

The  horizontal  grids  are  uniform.  The  grid  interval  is 
60  km.  The  pressure,  p,  is  used  as  vertical  coordinates  and 
k  is  the  vertical  grid  index.    The  pressure  of  Jrth     level, 


Pk 


in   the   model    is   specified   by   the   exponential 


relationship  with  height,  Zk , 


P  e 


zk 


(A.l) 


where  p0  is  the  pressure  at  the  model  bottom. 

Apply  (A.l)  to  the  total  height  of  the  model,  zt , 


:t  ■  "(-5T-)  • 


where  pt  is  the  pressure  at  model  top. 

Using   uniform   increment   of   height   in   vertical,    the 

corresponding  pressure  is 


(k-1) 


An(— 


■) 


po  e 


where  kt   is   the   total  number   of  vertical   levels  in   the 
model . 
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In  this  study,  the  total  number  of  vertical  levels,  kt , 
is  54.  The  corresponding  pressure  levels  used  in  the  model 
are  shown  in  Table  A.l 

TABLE  A.l  FORECAST  LEVELS  (cb) 


VEL 

p(ch) 

P(cb) 

P(cb) 

1 

95.86 

19 

44.81 

37 

20.95 

2 

91.9 

20 

42.96 

38 

20.08 

3 

88.1 

21 

41.18 

39 

19.25 

4 

84.45 

22 

39.48 

40 

18.25 

5 

80.96 

23 

37.85 

41 

17.69 

6 

77.61 

24 

36.28 

42 

16.96 

7 

74.4 

25 

34.78 

43 

16.26 

8 

71.32 

26 

33.34 

44 

15.59 

9 

68.37 

27 

31.96 

45 

14.94 

10 

65.54 

28 

30.64 

46 

14.32 

11 

62.83 

29 

29.37 

47 

13.73 

12 

60.23 

30 

28.16 

48 

13.16 

13 

57.54 

31 

26.99 

49 

12.62 

14 

55.35 

32 

25.88 

50 

12.10 

15 

53.06 

33 

24.81 

51 

11.60 

16 

50.87 

34 

23.78 

52 

11.12 

17 

48.76 

35 

22.80 

53 

10.66 

18 

46.75 

36 

21.85 

54 

10.22 
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B.   GRID  ARRANGEMENTS 

The  model  grids  and  the  all  field  variables  are  arranged 
as  in  Fig.  B.l.  The  vertical  level,  p,  is  the  midpoint 
between  levels  of  p.  Vertical  velocity  is  defined  at  p 
levels  and  the  other  variables  u,  v,  6,.  q  and  tt  are  defined 
at  p  levels.  The  horizontal  grids  are  defined  by  y  and  y. 
Variables  8,  q,  tt  and  (a  are  defined  at  y  where  co  j  ,  k  is 
defined  at  y  and  p  coordinates  and  is  calculated  from  the 
continuity  equation. 

"j,k  "    «j,k-l  -      .k      J"1       <vj>k-l  "  Vj-i,k-l) 

yj  -  yj-i 


yj-i      yj-i     yj        yj       yj+i     yj+i 


Pk+l 

Pk+l 
Pk 

Pk 
Pk-1 

Pk-1 


ej-l     luj-l 
---r-- 


T 


wJ,k 


ej+l       luj+l 


Fig.  B.l    The  model  grid  arrangement 
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C.   FINITE  DIFFERENCE  FORMULATION 

The  equation  set  (2.1.16)  and  ( 2 . 1 . 19 ) - ( 2 . 1 . 23  )  can  be 
integrated  by  a  pure  marching  process.  In  this  study  a 
spatially  staggered  finite-difference  scheme  is  used  as 
compared  to  the  unstaggered  scheme  used  in  Williams  (1974). 
The  variables  6,  ir  and  q  are  carried  at  the  same  points  in 
y,  and  u  and  v  are  carried  equidistant  between  these  points 
(see  Fig.  B.l).  Winninghoff  (1968),  Arakawa  and  Lamb  (1977) 
and  Schoenstadt  (1980)  have  shown  that  this  staggered  grid, 
which  is  known  as  scheme  B,  is  much  better  than  the 
unstaggered  grid  (scheme  A)  for  geostrophic  adjustment. 
This  scheme  is  especially  superior  when  small-scale  heating 
occurs,  such  as  near  a  front.  The  nonlinear  terms  are 
arranged  to  conserve  mean  squares  in  the  advection  terms. 
Leapfrog  time  differencing  is  used,  except  the  model  is 
restarted  every  fifth  timestep  with  an  Euler  backward  step 
following  Arakawa  and  Lamb  (1977). 

For  the  finite  difference  formulation,  we  discuss  only 
the  spatial  derivative  terms  in  the  u  and  v  equations  [i.e., 
(2. 1.16)  and  (2. 1. 19) ] , 


t-tv  vF]J,k   =   [<VJ+l,k  +  vJ,k)    (pj+l,k  +  Fj,k) 


ay 


-  (vj,k  +  vj-i,k)(Fj,k  +  Fj-i,kW(yj+i  -  yj> 
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C"KywF3J.*  =  CSj+lik+1  +  "j,k+l)    (Fj.ic+i  +  Fjjk) 

"   Bj+l,k  +  "j,k)    (Fj.k  +  Fjsk.1)]/4(pk+1  -  pk) 


tv^.k  =  (pj+i,k  -  Fj-i,k)/(yj+i  +  yj-i) 


2F  (FJ+l,k   "   Fj,k)  (Fjjk  -  Fj.1>k) 


•ty^j.k  =   C 


(yj+i  -  yj) 


/A  A 


(yj  -  yj-i) 


}/(yj+i  -  yj) 


_a       .8F_  r  Kk+1  (FJ,k+l  "  FJ.:0  MFj,k  -  Fj,k-On// 

Mb  K  kh,*  -  C : : r— ]/(Pk+l  -  Pk) 


(Pk+l  "  Pk) 


(Pk  "  Pk-l) 


where    F    is    either    u    or    v 
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For  the  9  and  q  equations, 


[-6T  vG]J»k  =  CvJ,k(GJ+l,k  +  Sj.k)  -  vj-l,k(Gj,k  +  Gj-i,k)/2(yj  -  yj-i) 


^Ty"6^.1*  =  ^j,k+l(Gj,k+l  +  Gj,k)   ""j,k(Gj»k  +  Gj,k-l)/2(Pk+l  -  Pk) 


frftl.k  =  fsj+i,k  -  Gj-i,k)/(yj+i  +  yj-i) 


^2G  (Gj+l,k  -  Gj,k) 

W*i*  =  L    (yi+i  -  y7) 


<Gj,k  "  Gj-l,k)     /A        . 

(yj  -yj-i)    ]/(yJ  *  yJ"15 


r5      „dG-,  -  Kk+1    (Gj,k+1   -  Gj,k) 

(Pk+1   "   Pk) 


MGj,k  "  Gj,k-l) 
(Pk  "  Pk-l) 


]/(Pk+l  "  Pk) 


where    G    is    either    0    or    q. 
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